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Executive  Summary 


The  principal  objective  of  tiie  present  research  proposal  v^^as  to  investigate, 
experimentally  and  numerically,  the  use  of  jets  to  actively  control  die  aerodynamic  forces 
for  high-speed  flight  vehicles.  As  a  test  bed,  single  port  injection  into  a  Mach  5.0 
crossflow  through  diamond  and  circular  orifices,  at  various  incidence  angles  and 
momentum  ratios  were  examined.  To  meet  the  objective,  a  research  program  was 
performed  to  characterize  and  control  the  secondary  flow  stmctures  associated  with  die 
jet  interaction  flows.  The  research  approach  was  divided  into  diree  main  thrusts.  First, 
fundamental  studies  to  characterize  die  mean  and  turbident  flow  structure  were 
performed.  Ihis  included  detailed  ©cperiments  at  Mach  5.0  using  advanced  optical 
diagnostics  and  numerical  simulations  using  two-equation  RANS  turbulence  modeling. 
The  flow  field  characterizaticm  included  shock  stmctures,  surface  pressure,  vorticity  and 
turbulence.  Second,  the  numerical  simulations  were  studied  to  assess  vorticity  and 
turbulent  transport  mechanisms.  Lastiy,  parametric  experiments  and  Reynolds  averaged 
Navier-Stokes  simidations  were  performed,  to  quantify  vorticity  control.  This  included 
manipulating  the  jet  shock  stmcture  to  modify  the  vorticity  field,  adding  downstream 
ramps  to  utilize  the  Magnus  force,  varying  the  injector  geometry,  incidence  angle  and 
momentum  ratio. 

The  Air  Force  Research  Laboratory  High-Speed  Propidsion  Division  is 
investigating  the  use  of  diamond  jets  for  supersonic  combustion  ramjet  (scramjet) 
primary  fuel  injectors. 


The  results  of  this  research  program  include  a  new  experimental  and  numerical 
database  for  sonic  injection  into  a  Mach  5.0  crossflow  and  new  understanding  of  the 
vorticity  dynamics  for  this  flow.  The  results  were  disseminated  to  the  US  Air  Force  and 
aerospace  industry  through  this  contractor  report,  archival  journal  publications  and 
presentations  at  national  conferences. 

The  lead  PI  for  this  project  was  Dr.  Rodn^  D.  W.  Bowersox,  Associate  Professor 
of  Aerospace  Engineering  at  Texas  A&M  University.  He  moved  to  Texas  A&M 
University  during  the  second  year  of  the  project,  and  the  project  remained  at  The 
University  of  Alabama  under  the  co-PI  Dr.  Charles  L.  Karr.  Dr.  Muhammad  A.  R.  Sharif 
was  die  third  co-PI  on  diis  project.  The  Mach  5.0  wind  tunnel  tests  and  numerical 
simulations  were  carried  out  by  Dr.  Bowersox  and  his  graduate  students  Huaiguo  Fan 
(PhD,  currently  on  the  faculty  at  the  Chongqing  University),  Daina  Lee  (MS,  currently  at 
ATA,  Inc.  at  AEDC)  and  Ravichandra  Srinivasan  (currently  a  PhD  candidate  at  Texas 
A&M).  Drs.  Karr  and  Sharif  performed  the  design  optimization  studies.  Althou^, 
unfunded.  Dr.  Bowersox  continued  die  numerical  simulations  at  Texas  A&M  University. 
In  addition,  he  is  currendy  performing  detailed  high-fidelity  simulation  under  ai  Air 
Force  supercomputer  grant  on  die  Air  Force  ASC  HPC;  the  results  from  these  simulations 
will  be  reported  at  a  later  date. 
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Chapter  1 

Introduction 


1.1  Documented  Air  Force  Requirements 

Tlie  U.S.  Air  Force  requirements  for  hypersonic  vehicles,  and  the  tedmical  diallenges 
diereof,  are  cleariy  outlined  in  (1)  die  New  World  Vista’s  document  (1995),  (2)  the 
System  Concepts  and  Technology  Development  of  Operationally  Responsive  Reusable 
Military  Aerospace  Vehicle  p^r  [Dolvin  (2002)],  and  (3)  the  Unsteady  Aerodyrumics 
and  Hypersonics  section  of  the  AFOSR  BAA  [Schmisseur  (2003)].  In  summary, 
“sustained  hypersonic  fli^t  offers  potential  revolutionary  improvements  in  future  war- 
fighting  and  space  laundi  c^abilities.  ...  Hypersonic  speeds  greatly  enhance  bodi 
survivability  and  lethality  by  practically  eliminating  enemy  defensive  capabilities.  (Vew 
World  Vistas,  1995).”  Potential  future  Air  Force  vehicles  include  space  launch  support 
accelerators,  missiles,  maneuvering  re-entry  vehicles,  and  r^id-response  global-reach 
aircraft  systems. 


Fig.  1.1  Example  Air  Force  Hypersonic  Vehicles  [Erbland  (2002)] 
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Many  of  flie  advanced  hypersonic  vehicle  concepts,  such  as  those  in  Fig.  1.1,  are 
envisioned  to  use  complex  combined  engines  and/or  perform  rapid  maneuvers,  such  as 
that  shown  in  Fig.  1.2 


Fig.  1.2  Generic  Scramjet  Engine 


High-speed  jet  interaction  flows,  as  shown  in  Fig.  1.3,  have  numerous  potential 
applications  for  die  vdiicles  shown  in  Fig.  1.1.  For  trample  jets  can  be  very  effective 
rapid  reaction  control  devices  for  hypersonic  vehicles.  One  advantage  of  a  reaction 
control  jet  over  a  ctMitrol  sur&ce  is  the  elimination  of  a  thermally  protected  control 
surface.  However,  jet  interaction  flows  are  turbulent  (see  Fig.  1.3a)  and  di^  are 
characterized  by  numerous  complicated  vortical  flow  features,  as  depicted  in  Fig.  1.3b. 
Thus,  nodem  physical  models  are  not  accurate.  Jet  flows  and  prediction  limitation  are 
described  in  more  detail  in  Chapter  2. 
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(a)  Mie  Scattering  Image  (b)  Numerical  Simulation 

Fig.  1.3  Sonic  Injection  into  a  Mach  5  Cross  Flow  [Present  Study] 

Jets  are  also  effective  for  film  cooling  fiiermal  protection.  A  Mach  8  vehicle  will 
experience  wall  temperatures  on  tiie  order  of  2,700K.  The  only  metals  that  can  survive 
these  tenqieratures  are  very  heavy  (e.g.,  Tungstei  is  nearly  ten  times  heavier  than 
Titairium).  Thus,  the  aerosurfaces  will  require  thermal  protection  (either  ceramics  or  film 
cooling).  One  of  the  distinguishing  characteristics  of  a  jet  into  a  cross- flow  is  the  strong 
counter-rotating  vortex  pair  within  the  plume  (see  Chaptw  2  for  a  detailed  flow  field 
description).  This  counter  rotating  vortex  pair  has  an  adverse  influence  on  the  film 
cooling  effectiveness,  where  the  vortex  pair  interacts  with  the  surface  causing  the  jet 
plume  to  lift-off  tire  surface  via  the  Magnus  force.  In  a  recent  stiufy  aimed  at  improving 
the  paietration  of  a  jet  into  a  supersoruc  cross- flow,  our  group  investigated  the  use  of 
small  ramps  placed  just  downstream  of  the  injector  port  [Wilson,  Bowersox  and  CHawe 
(1998)].  The  ramps  were  designed  to  interact  with  the  counter-rotating  vortex  pair,  and 
thus  uiCTease  the  penetratirm  via  the  Magnus  force.  Our  team  discovered  that  by  varying 
tire  ramp  geometry  the  plume  vortidty  could  be  controlled.  One  ramp  configuration 
reduced  the  vorticity  within  the  plume  to  the  level  of  being  imdetectable,  and  our 


numerical  simulations  [Chenault,  Beran  and  Bowersox  (1999)]  indicated  that  the 
recompression  shock  (described  in  Chapter  2)  is  a  main  vorticity  production  mechanism. 
Hence,  the  most  likely  case  for  the  plume  vorticity  control  is  that  the  ramps  reshape  the 
recompression  shock.  Once  the  flow  processes  are  better  understood,  mefliodologies  that 
range  from  simple  fixed  ramps  to  sophisticated  nonlinear  adaptive  control  using  micro- 
electrical- mechanical  systems  (MEMS)  can  be  developed  to  control  the  jet  plume 
features  including  plume  axial  vorticity. 

Although  the  counter-rotating  vortex  pair  can  have  a  negative  impact  on  film 
cooling,  this  characteristic  has  the  potential  to  be  an  effective  means  of  boundary  l^er 
control.  It  is  expected  that  the  axial  vorticity  widiin  the  jet  plume  can  be  used  to  energize 
the  lower  portion  of  the  boundary  layer  wifii  the  goal  of  reducing  wall  pressure 
fluctuations  by  extending  the  axial  length  of  the  "U- vortex"  described  in  Hinze  (1975). 
Low-speed  experimental  data  indicate  that  the  wall  shears  stress  is  linearly  proportional 
to  the  root-mean- square  wall  pressure  fluctuations  (Hinze,  1975).  Hence,  the  reduced 
wall  pressure  fluctuations  will  result  in  reduced  momentum  loss  (skin  fiiction). 
Energizing  the  boundary  layer  can  also  be  used  to  increase  the  resistance  of  the  boundary 
layer  to  shock  induced  separation.  These  include  reaction  control,  separation  control, 
shock  interaction  control  and  fuel  injection  for  the  scramjet  engine  as  shown  in  Fig.  1.2. 

In  summary,  the  understanding  and  predictability  of  a  jet  into  a  cross-flow, 
especially  at  high- speeds,  is  limited  by  the  extreme  dearth  of  experimental  studies 
(Chapter  2  describes  in  detail  the  current  limitations).  However,  as  discussed  Chapter  2, 
recent  advances  in  experimental  methods  allow  for  in  depth  studies  directed  at 
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imderstanding  the  flow  physics.  The  better  flow  field  understanding,  produced  in  this 
stu(fy,  provides  engineers  and  scientists  the  requisite  information  to  capitalize  on  the  jet 
flow  features  to  enhance  die  performance  of  hypersonic  systems. 

13  Research  Objectives  and  Approach 

The  principal  objective  of  the  present  research  proposal  is  to  investigate, 
experimentally  and  numerically,  the  use  of  jets  to  actively  control  the  aerodynamic  forces 
for  high-speed  flight  vehicles.  As  a  test  bed,  single  port  injection  into  a  Mach  5.0 
crossflow  through  diamond  and  circular  orifices,  at  various  incidence  angles  and 
momaitum  ratios  were  examined.  To  meet  the  objective,  a  research  program  was 
performed  to  characterize  and  control  the  secondary  flow  stmctures  associated  with  the 
jet  interaction  flows.  The  research  approach  was  divided  into  three  main  thrusts.  First, 
fundamental  studies  to  characterize  the  mean  and  turbulent  flow  structure  were 
performed.  This  included  detailed  experiments  at  Mach  5.0  using  advanced  optical 
diagnostics  and  numerical  simulations  using  two-equation  RANS  turbulence  modeling. 
The  flow  field  characterization  included  shock  structures,  surface  pressure,  vorticity  and 
turbulence.  Second,  the  numerical  simulations  were  studied  to  assess  vorticity  and 
turbulent  transport  mechanisms.  Lastly,  parametric  experiments  and  Reynolds  averaged 
Navier-Stokes  simulations  were  performed  to  quantify  vorticify  control.  This  included 
manipidating  the  jet  shock  structure  to  modify  the  vorticify  field,  adding  downstream 
ramps  to  utilize  the  Magnus  force,  varying  the  injector  geometry,  incidence  angle  and 
momentum  ratio. 
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1.4  Scientific  Impact  and  Technology  Transitions 


The  Air  Force  Research  Laboratory  High-Speed  Propulsion  Division  is 
investigating  the  use  of  diamond  jets  for  supersonic  combustion  ramjet  (scramjet) 
primaiy  fuel  injectors. 

1.5  Deliverables 

The  results  of  this  research  program  include  a  new  experimental  and  numerical 
database  for  sonic  injection  into  a  Mach  5.0  crossflow  and  new  imderstanding  of  tiie 
vorticity  ^namics  for  this  flow.  The  results  were  disseminated  to  the  US  Air  Force  and 
aerospace  industry  through  this  contractor  rqrort,  archival  journal  publications  and 
presentations  at  national  conferaices. 

1.6  Personnel 

The  lead  PI  for  this  project  was  Dr.  Rodney  D.  W.  Bowersox,  Associate  Professor 
of  Aerospace  Engineering  at  Texas  A&M  University.  He  moved  to  Texas  A&M 
University  during  the  second  year  of  die  project,  and  the  project  remained  at  The 
University  of  Alabama  under  the  co-PI  Dr.  Charles  L.  Karr.  Dr.  Muhammad  A.  R.  Sharif 
was  die  third  co-PI  on  this  project.  The  Mach  5.0  wind  tunnel  tests  and  numerical 
simulations  were  carried  out  by  Dr.  Bowersox  and  his  graduate  students  Huaiguo  Fan 
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(PhD,  currently  on  the  faculty  at  the  Chongqing  University),  Daina  Lee  (MS,  currently  at 
ATA,  Inc.  at  AEDC)  and  Ravichandra  Srinivasan  (currently  a  PhD  candidate  at  Texas 
A&M).  Drs.  Karr  and  Sharif  performed  tiie  design  optimization  studies. 

1.7  Publication  Summary  (*indicates  student  under  die  direction  of  Dr.  Bowersox) 

Bowersox,  R.,  Fan,*  H.,  and  Lee,*  D.,  “Sonic  Injection  into  a  Mach  5.0  Freestream 
throu^  Diamond  Orifices,”  Journal  cf  Propulsion  and  Power,  Vol.  20,  No.  2,  2004,  pp. 
280-287. 

Srinivasan,*  R.  and  Bowersox,  R  “Use  of  Transverse  Injection  for  Hypersonic  Flow 
Control,”  AIAA-2004-2698, 2nd  AIAA  Flow  Control  Conference,  Portland,  Oregon  Jun  28 
-  Jul  1, 2004. 

Bowersox,  R.,  Fan*,  H.,  and  Lee*,  D.,  “Sonic  Injection  into  a  Mach  5.0  Freestream 
dirough  Diamond  Shaped  Orifices  at  Various  Incidence,”  AIAA-2003-5190, 

AIAA/ASME/SAE/ASEE  Joint  Propulsion  Conference,  Huntsville,  AL,  July  2003. 

Srinivasan,  R.  and  Bowersox,  R.,  “Numerical  Analyses  of  Gaseous  Injection  into  a 
Mach  5.0  Freestream  through  Diamond  Orifices,”  AIAA-2003-5188, 

AIAA/ASME/SAE/ASEE  Joint  Propulsion  Conference,  Huntsville,  AL,  July  2003. 

Fan*,  H.  and  Bowersox,  "Gaseous  Injection  through  Diamond  Orifices  at  Various 
Incidence  Angles  into  a  Hypersonic  Freestream,"  AIAA  2001-1050,  39^**  AIAA 
Aerospace  Sciences  Meeting,  Reno  NV,  Jan  2001. 
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Chapter  2 

Jet-Interaction  Background  Review 
2.1  Jet-Into-A-Cross-flow 

Injection  into  a  cross-flow  is  a  fundamentally  important  flow  field  diat  has 
implications  across  a  myriad  of  military  and  civil  disciplines.  Thus,  numerous 
experimental  studies,  primarily  at  low-speeds,  have  been  performed  to  quantitatively 
assess  the  mean  flow  field  and  qualitatively  examine  the  instantaneous  vortical  flow 
features  [for  example,  Margason  (1993),  Kuethe  (1935),  and  Kamotani  and  Greber 
(1972)].  However,  studies  that  contain  turbulence  data  suitable  for  flow  field 
understanding  are  exceedingly  scarce  [Margason  (1993)].  Furfliermore,  the  available 
mean  and  turbulence  data  for  supersonic  flows  are  also  limited.  However,  some  studies 
were  found,  and  the  results  of  which  are  integrated  into  die  following  the  discussions. 
Hypersonic  jet- interaction  studies  are  virtually  non-existent. 

Injection  into  a  low-speed  cross-flow  produces  a  complex  flow  with  at  least  four 
distinguishable  vortical  flow  characteristics  (see  Fig.  2.1).  As  the  jet  emerges  into  the 
cross-flow,  it  is  curved  downstream  by  the  cross-flow.  The  jet  plume  cross-section 
evolves  into  a  cardioidal  or  “kidney-bean”  shape.  The  counter-rotating  vortex  pair  within 
the  plume  is  responsible  for  the  distorted  plume  shape.  A  number  of  flow  mechanisms 
are  responsible  for  generating  this  secondary  flow.  For  example,  turning  of  the  flow  into 
die  downstream  direction  creates  a  pressure  gradient  across  the  flow  that  induces  vortical 
motion.  Lateral  shearing  along  the  plume  edges  also  contributes  to  the  vorticity.  The  jet- 
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free  stream  interaction  (separation)  creates  a  horseshoe  vortex  similar  to  that  of  a  wing- 
body  junction  on  an  aircraft.  The  third  vortex  system  consists  of  jet-shear  I^er  vortices. 
Lastly,  a  fourth  wake  vortex  system  exists.  Until  recently,  the  formation  of  this  unsteady 
vortex  phenomaion  had  been  attributed  to  vortex  shedding  similar  the  Karman-vortex 
wake  of  a  cylinder  in  a  cross-flow.  Because  Fric  and  Roshko  (1994)  were  only  able  to 
visualize  this  system  when  the  boundary  layer  was  seeded  with  smoke,  th^  reasoned  that 
the  vortices  were  the  result  of  upstream  boundary  layer  vorticity  being  transported  up  into 
the  wake.  The  process  is  believed  to  be  the  result  of  “separation  events”  and  eruptions  of 
the  boundary  layer  fluid  and  vorticity  into  “tornado- like”  structures  that  begin  at  the  wall 
and  are  entrained  into  the  jet  flow.  Although  Fric  and  Roshko  present  rather  convincing 
evidence  for  this  hypothesis,  some  contradicting  evidence  has  also  been  presented.  For 
example,  the  planar  laser  induced  fluorescence  concentration  measurements  of  Lozano  et 
al.  (1994),  show  signs  of  the  structures.  However  in  that  sturfy,  only  the  jet  flow  was 
seeded.  Hence,  these  data  may  suggest  that  these  structures  could  in  part  emanate  from 
the  jet  and  not  the  wall.  Adding  to  this  apparent  dilemma,  Kelso  et  al.  (1993)  present 
three  possible  scenarios  leading  to  the  formation  of  the  wake  vortex  system  These 
unsteaxfy  phenomena  can  be  very  important  for  future  flow  control  algorithms. 
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For  high-speed  flows,  the  mean  flow  features  are  similar  to  those  described  above. 
However,  compressibility  creates  additional  features  that  do  not  have  incompressible 
counterparts.  For  example,  perpendicular  injection  of  an  underexpanded  sonic  or 
supersonic  jet  into  a  si5)ersonic  free  stream  produces  several  flow  structures.  The  first  of 
these  is  a  bow  shock  produced  as  the  free  stream  impacts  the  injection  streamtube;  in  fiiis 
respect  the  injectant  acts  as  a  solid  cylindrical  bo<fy  [Hollo,  McDaniel  and  Hartfield 
(1994)].  For  injector  configurations  where  injector  diameter  is  greater  dian  the  boundary 
layer  diickness,  a  separation  bubble  and  a  lambda  shock  form  slightly  upstream  of  die 
injector  port  [Schetz  and  Billig  (1966)].  After  entering  the  free  stream,  the  jet 
experiences  a  rapid  Prandti-Meyer  expansion  (usually  assumed  to  be  an  isentropic 
process)  surrounded  by  a  barrel  shock  [Grasso  and  Magj  (1995)].  A  shock  normal  to  the 
jet  path,  known  as  die  Mach  disk,  terminates  the  barrel  shock,  and  compresses  the  flow  to 
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tile  appropriate  (i.e.,  effective)  back  pressure.  As  the  free  stream  flow  wraps  around  the 
jet,  a  recompression  shock  is  created.  Recent  CFD  calcidations  by  our  team  [Chenault, 
Beran  and  Bowersox  (1999)]  (see  Fig.  2.2)  indicate  that  the  recompression  shodk:  creates 
an  addition  counter-rotating  vortex  pair  that  eventually  engulfs  the  original  jet  induced 
pair  discussed  above. 


Fig.  2.2  High-Speed  Jet-Interaction  Vortidty  Distribution  [Chenault,  Beran  and  Bowersox  (1999)] 


An^ed  injection  is  a  means  of  reducing  total  pressure  loss.  In  either  normal  or 
angled  injection,  the  entry  of  the  injectant  jet  into  the  mainstream  flow  can  be  regarded  as 
a  two-stage  process  [Schetz  and  Billig  (1966)].  The  jet  first  enters  the  main  flow  and 
remains  relatively  intact  as  it  sqiands  to  tiie  height  of  the  Mach  disk.  Beyond  the  Mach 
disk,  the  flow  turns  coaxial  to,  and  accelerates  witii,  tiie  main  flow.  In  the  second  stage, 
the  jet  acts  as  a  coaxial  vortex  mixing  structure.  It  is  this  feature  that  may  be  useful  for 
boundary  layer  control  (i.e.,  drag  reduction). 

During  a  rather  extensive  literature  review,  only  a  very  limited  number  of  studies 
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on  injection  flows  were  found  that  included  turbulence  measurements.  For  low-speed 
flows  Andreopoulus  and  Rodi  (1984),  Fric  and  Roshko  (1989),  Kamotani  and  Greber 
(1972)  were  the  only  ones  identified.  For  high-speed  flows,  the  McCann  and  Bowersox 
(1996)  and  Bowersox  (1996c)  Mach  3.0  studies  were  the  only  ones  located.  This  finding 
is  consistent  with  that  of  the  Margason  (1993).  Hence,  the  turbulent  flow  physics  of  the 
injection  into  a  cross-flow  are  significantly  less  understood  than  their  mean  flow 
counterpart  However  from  file  available  data,  some  limited  qualitative  insight  into  the 
overall  structure  can  be  gleaned.  For  example  the  axial  turbulence  intensity  contour  plots 
of  Kamotani  and  Greber  (1972)  indicated  that  two  peaks  in  the  intensity  levels  existed 
across  the  plume.  In  addition,  the  peaks  were  roughly  co- located  with  each  of  the 
vortices  of  the  counter  rotating  vortex  pair.  Investigating  a  high-speed  angled  injection 
case,  McCann  and  Bowersox  (1996)  found  fiiat  the  turbulent  kinetic  energy  also  had  two 
peaks,  which  were  roughly,  located  just  below  each  vortex  with  the  vortex  pair.  They 
deduced  tiiat  the  increased  strain  rates  due  to  fiie  secondary  vortex  flow  resulted  in  an 
increased  production  of  the  turbulent  kinetic  energy.  The  data  of  McCann  and  Bowersox 
also  demonstrated  that  the  x-y  turbulent  shear  stresses  peaked  along  the  plume  centerline, 
and  fire  x-z  shear  stresses  were  highly  asymmetrical.  However,  as  suggested  by  Cantwell 
(General  Discussion,  1993),  the  peak  turbulence  levels  along  the  jet  centerline  may  be  the 
result  of  an  unsteady  movement  of  the  vortex  pair.  It  is  also  important  to  mention  here 
that  a  number  of  experimental  studies  show  an  asymmetiy  in  the  mean  flow.  Further,  the 
numerical  simulation  of  a  shocked  jet  (Drummond,  1996)  also  showed  signs  of  an 
asymmetry.  Bowersox  (1996c)  presents  three  possible  scenarios  to  explain  the 
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generation  of  the  asymmetry.  The  incompressible  data  of  Andreopoulus  and  Rodi  (1984) 
presented  flie  axial  variation  of  the  x-y  and  xz  turbulent  shear  stresses  for  two  spanwise 
locations.  Unfortunately,  the  spatial  resolution  of  that  study  was  sparse;  hence,  only 
limited  quantitative  information  was  available.  However,  those  data  showed  trends  in  the 
turbulence  dec^  rates.  Though  the  limited  data  available  in  the  literature,  in  conjunction 
with  a  current  understanding  of  the  mean  flow  field,  can  provide  a  qualitative  impression 
of  the  turbulent  flow  structure,  studies  with  enough  detail  for  turbulence  model 
development  and  validation  are  virtually  non-existent.  This  conclusion  is  supported  by 
the  review  of  Margason  (1993),  and  it  was  also  the  consensus  during  the  General 
Discussion  session  following  a  recent  AGARD  meeting  (General  Discussion,  1993) 
among  a  number  of  key  international  researchers.  Further,  Chui  et.  al.  (1993)  and 
Alverez  et  al.  (1993)  both  indicated  that  turbulence  modeling  was  the  limiting  key  factor 
in  determining  the  accuraty  of  their  numerical  predictions.  Hence,  detailed  experimental 
studies  are  required  to  develop  turbidence  models  for  design  applications. 

2.2  Overview  of  Current  Prediction  Methods  and  Empirical  Database 

Currently,  engineers  and  scientists  rely  on  an  approximate  form  the  governing 
equations  of  motion  to  predict  the  overall  mean  flow  character  of  these  flows.  The 
Reynolds  and  Favre  averaged  forms  of  the  Navier-Stokes  equations  are  universally  used 
to  compute  low^speed  and  high-speed  flows,  respectively. 

Associated  with  the  averaging  procedure  is  the  requirement  for  accurate 
turbulence  models.  Currently,  industry  standard  eddy  viscosity  gradient  transport 
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turbulence  models,  such  as  Baldwin-Lomax  algebraic  and  k-s/co  two  equation,  are 
extended  to  more  complex  flows  on  an  orf  hoc  basis.  However,  as  discussed  in  Wilcox 
(1993),  current  eddy  viscosity  models  are  not  accurate  for  three-dimensional,  vortex- 
dominated  flows.  As  discussed  in  Section  2.1,  the  jet  into  a  cross- flow  flow  field  is 
dominated  by  numerous  three-dimensional  vortex  systems;  thus,  the  current  prediction 
methods  are  not  accurate. 

In  a  recent  numerical  study  of  injection  into  a  cross-flow,  Chui  et  al.  (1993)  noted 
that  the  predictions  of  the  flow  with  the  eddy  viscosity  Baldwin-Lomax  and  Baldwin- 
Barth  models  "...  compared  no  better  with  the  experimental  data  than  the  laminar  flow 
computation.”  Alverez  et  al.  (1993)  reported  similar  shortcomings.  In  a  recent  survey  of 
333  articles  that  covered  a  50-year  era,  Margason  (1993)  concluded  fliat  there  is  a  need 
for  high  quality,  high  fidelity  experimental  data  that  will  allow  for  verification  of  current 
and  future  computational  fluid  dynamic  results  and  to  define  the  unsteady  flow 
characteristics.  Also,  the  during  the  General  Discussion,  following  the  1993  AGARD 
(General  Discussion,  1993)  conference,  entitled  Computational  and  Experimental 
Assessment  of  Jets  in  a  Cross-flow,  diere  was  some  surprise  that  more  detailed  and  high 
fidelity  Laser  Doppler  and  Particle  Image  Velocimetiy  data  were  not  available  for 
turbulence  model  development  and  evaluation. 

2.3  Recent  Advances  in  Technology  that  Allow  for  the  Present  Research 

Modem  advanced  laser  based  experimental  diagnostic  techniques,  such  as  particle 
image  velocimetiy  and  pressure  sensitive  paint,  allow  for  detailed  and  accurate 
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measurements  of  the  instantaneous  mean,  and  turbulence  statistics  at  high  spatial 
resolution.  In  addition,  modem  computational  resources,  such  as  the  University  of 
Alabama  Super  Computer  (described  in  Chapter  3),  allow  for  routine  simulation  of 
complex  three-dimensional  flows. 

2.4  Control  of  the  Jet  into  a  Cross-flow  Using  Variable  Wall  Geometry 

In  a  recent  research  program  focused  on  supersonic  combustion  applications, 
Wilson,  Bowersox  and  Glawe  (1998)  placed  small  ramps  just  downstream  of  the  exit  port 
of  a  jet  issuing  into  a  high-speed  cross- flow.  The  original  motivation  behind  placing  the 
ramps  was  to  increase  penetration  via  tire  Magnus  force  garerated  by  interaction  between 
the  strong  counter  rotating  vortex  pair  within  the  plume  (Section  2.1  describes  the  flow 
features)  and  the  ramp  surfaces.  That  study  was  limited  to  one  injector  angle  (25 
degrees)  at  Mach  3.0.  However,  seven  ramp  configurations  were  tested.  Because  of  the 
lack  of  design  tools  (i.e.,  suitable  turbulence  models),  the  ramps  were  designed  based  on 
intuitive  arguments.  The  effects  of  the  small  ramps  on  the  plume  structure  were 
significantly  more  pronounced  than  expected,  and  they  were  strongly  dependent  on  the 
ramp  geometry.  For  one  particular  amp,  the  normally  very  strong  counter  rotating 
vortex  pair  (20,000-30,000  /s)  was  reduced  in  strength  to  the  point  of  being  undetectable 
with  laser  based  Mie  scattering  experimental  apparatus.  On  the  other  hand,  a  second 
ramp  configuration  produced  stronger  vortices  within  the  plume.  The  aforementioned 
lack  of  understanding  of  the  flow  physics  made  determining  the  physical  mechanisms 
responsible  for  the  effects  of  the  ramps  on  the  plume  structure  unattainable.  However, 
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the  recent  simulations  of  Chenault,  Beran  and  Bowersox  (1999)  indicated  that  the 
recompression  shock  is  a  main  vorticity  production  mechanism.  Hence,  the  most  likely 
case  for  the  plume  vorticily  control  is  that  the  ramps  reshape  the  recompression  shock. 
Once  these  flow  processes  are  better  understood,  methodologies  diat  range  from  simple 
fixed  ramps  to  sophisticated  nonlinear  adaptive  control  using  micro-electrical- mechanical 
systems  (MEMS)  can  be  developed  to  control  the  jet  plume  features  (e.g.,  shape, 
trajectory  and  unsteady,  three-dimensional  nonlinearities).  Thus,  optimizing  the 
performance  of  numerous  aigineering  endeavors. 
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Chapter  3 

Facilities  and  Apparatus  Development 

This  section  provides  a  detailed  description  of  die  facility  and  apparatus 
development  woA  diat  was  accomplished  in  the  research  {Hogram. 

3.1  Wind  Tunnels 

3.1.1  Mach  5. 0  Wind  Tunnel 

A  Mach  5.0  blow-down  tmmel  with  a  6.0  second  run  time  was  used  for  the 
present  study.  Hie  ftee  stream  Mach  number  was  measured  as  4.9  with  a  +2%  variation 
across  the  test  section  core  flow.  The  operating  conditions  are  listed  in  Table  3.1.  A 
schematic  of  die  wind  tunnel,  along  widi  a  photograph,  is  given  in  Fig.  3.1. 


Table  3.1  Freestream  Flow  Conditions 


Mach 

Pi(MPa) 

2V(ig 

Re/m  (xlO*) 

S  (mm) 

4.9512.0%* 

2.4044.0%  * 

360±3.0%‘ 

53. 

7.6±7.0%* 

“Variation  across  the  freestream,  constant  to  within  this  amoimt, 
^Variation  across  the  tunnel  span 
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(a)  Facility  Schematic 


(b)  Photographs 
3.1  Mach  5.0  Wind  Tunnel 


Hie  high-pressure  air  supply  pumping  system  was  a  Hypres  Air  System  Model 
HP-6000-NA4-E3  four-stage,  reciprocating  air-cooled  air  compressor,  which  was  driven 
by  a  7.46  kW  (10  hp),  230  VAC  Marathon  Electric  motor.  This  system  provides  a  flow 
rate  of  5.66  liter/s  (12  cfin).  The  air  from  the  compressor  was  cooled  and  dried  by 
inerstage  coolers  and  separators,  and  purifier  desiccant  cartridges.  The  system  provided 
grade  E  breathing  air  [CO  <  10  ppm  v/v,  CO2  <  500  ppm  v/v,  total  Itydrocarbon  content  < 
25  ppm  v/v,  oil  (condensable  mist/vi^or)  <  5  ppm  w/v]  with  a  dew  point  of  -45  °C  [63 
ppm  v/v  (0.05  mg/1)].  Ihe  30  MPa  (4400  psia)  high-pressure  air  was  stored  in  eight  DOT 
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air  bottles,  each  with  a  volume  of  0.058  n?.  The  time  to  compress  from  atmospheric  to 
full  capacity  was  approximately  four  hours,  and  a  fiill  supply  of  air  was  sufficient  for 
nominally  four  wind  tunnel  runs.  A  second  Hypres  Brand  H*-6000-NA3-E3  Ihree-stage 
compressor  system  was  purchased.  This  Ihree-stage  system  is  very  similar  to  the  above- 
described  four-stage  system.  The  system,  which  driven  by  a  5.59  kW  (7.5  hp)  motor 
provides  an  additional  4.25  liter/sec  (9  cfin)  grade  E  air. 

Because  the  mass  flow  rate  for  the  tunnel  exceeded  that  available  directly  from 
the  bottle  farm,  the  air  for  single  run  was  transferred  to  an  intermediate  tank  for 
temporary  storage  and  high  volume  blow-down.  The  intermediate  tank  consisted  of  three 
6.0  m  sections  of  10.16  cm  (4  inch)  diameter  schedule  120  stainless  pipe.  The  tank 
volume  was  0.11  nJ,  and  the  intermediate  pressure  was  set  at  15.9  MPa  (2,300  psia). 
The  air  from  the  intermediate  tank  flowed  to  the  tunnel  tiirough  a  2.54  cm  (1.0  inch) 
stainless  steel  tube.  Two  Grove  Model  302G  Pressure  Reducing  regulators,  located  in 
series,  regulated  the  pressure  of  the  air  entering  the  wind  tunnel  to  within  ±2%  of  the 
desired  pressure.  The  dome  pressure  for  the  first  regulator  was  set  at  5.86  MPa  (850  psia) 
and  the  second  was  set  5.17  MPa  (750  psi).  These  settings  overcame  the  pressure  drop  in 
the  system  (plumbing  and  pebble  bed  heater)  to  provide  a  total  pressure  at  the  nozzle  of 
2.34  MPa  (340  psi). 

To  prevent  oxygen  liquefaction,  a  pebble  bed  was  designed  to  heat  the  air  before 
entering  the  nozzle.  The  pebble  bed  heating  system  was  designed  and  constructed  in- 
house  for  this  project.  A  rectangular  (interior  dimensions  of  12.7  x  12.7  x  49.53  crrf) 
high  temperature  grade  stainless  steel  (RA321)  vessel  held  the  pebble  bed.  The  pressure 
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vessel  was  constructed  from  5.08  cm  diick  plates,  and  was  pressure  tested  widi  hydraulic 
oil  to  17.24  MPa  (2500  psi).  Shown  in  Fig.  3.2  is  a  photograph  of  die  pebble  bed.  The 
pebble  bed  consisted  of  stainless  steel  (316)  pebbles  with  diameter  of  6.35mm  each.  The 
pebbles  absorbed  thermal  energy  when  low-pressure  heated  air  flowed  through. 


The  pebbles  were  heated  widi  two  Reheat  Brand  model  PFO-12  all  stainless 
electric  air  heaters  (1.2  kW  each).  An  Eagje  Model  C7180V1  reciprocating  single  stage 
compressor  [5.4  kW  (7.5  hp)]  provided  die  low-pressure  (0.2  MPa)  air  at  a  rate  of  9.1 
n^/min  (30  cfrn).  An  Arrow  brand  model  RH203  headess  twin  tower  regenerative 
desiccant  air  dryer  (40°  C),  with  micron  particulate  filters,  was  used  to  ensure  diat 
moisture  and  oil  v^or  did  not  accumulate  on  the  pebbles  during  preheat.  The  pr^at 
system  required  two-and-a-half  hours  to  heat  the  pebble  bed  fiem  room  temperature  to 
367  K  (200°F).  The  higji-pressure  air  for  a  tunnel  run  obtained  thermal  energy  via 
convection  from  the  pebbles  when  passing  through  the  pebble  bed. 
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The  Mach  5  contoured  nozzle  was  34.3  cm  long  with  Ihroat  area  of  1.93  cnf.  A 
photograph  of  the  nozzle  is  presented  in  Fig.  3.3.  The  nozzle  was  designed  by  the  PI  and 
constructed  in  1995  from  a  high  temperature  grade  (type  347)  stainless  steel.  The 
original  flow  straightner  consisted  of  a  single  aerogrid  (i.e.,  porous  plate)  with  forty  6.35 
mm  diameter  holes  located  approximately  3.8  cm  upstream  of  the  nozzle  contraction. 
The  Mach  number  was  uniform  to  within  2.0%  across  the  free  stream.  However,  the 
Mach  number  in  the  boimdary  layer  varied  significantly  in  the  span  wise  direction. 
Because  of  this  variation,  a  new  flow  straightner  was  designed  and  constructed.  The  new 
straightner  consisted  of  two  aerogrids  each  with  sixty  3.2  mm  holes  located  2.54  cm 
upstream  of  a  new  turbulence-dampening  screen.  Figure  3.4a  shows  the  new  flow 
straightner  and  the  original  aerogrid.  The  new  straightner  was  constructed  from  a  high 
temperature  grade  (type  RA321)  stainless  steel.  The  screen  was  located  approximately 
5.1  cm  upstream  of  the  nozzle  contraction.  It  is  also  noteworthy  that  the  stilling- 
chamber-to-throat  contraction  area  ratio  for  increased  from  13;1  to  19.4:1  with  the  new 
strai^tner.  Presented  in  Figs.  3.4b  and  c  is  a  comparison  of  the  flow  before  and  after  the 
new  flow  straightner  measured  with  the  Pitot  tube  in  the  five-hole-probe.  The  boundary 
layer  with  the  new  straightner  is  significantly  more  uniform  across  the  span  of  the  flow. 
However,  the  free  stream  was  only  modestly  improved. 
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(a)  Sidewall  Removed  with  Old  Aerogrid 


Fig  3.3  Photographs  of  the  Mach  5.0  Contoured  Nozzle 

The  stilling  chamber  Pitot  pressure  was  measured  widi  a  0-3.4  MPa  gage  Endevco 
Brand  Model  8510B-500  piezoresistive  pressure  transducer  widi  sensitivity  8.44 
mvolt/pa  The  miniature  pressure  transducer  (3.8  mm  face  diameter)  had  a  resonance 
frequency  of  500  kHz  and  stable  operation  over  a  temperature  range  of  0-366  K.  It  was 
connected  to  a  2.38  mm  outer  diameter  stainless  tube,  which  was  inserted  into  the 
chamber  and  mounted  to  the  nozzle  sidewall  by  a  high-temperature  resistant  epoxy  seal. 


The  transducer  was  connected  to  an  Bidevco  Model  4428A  Pressure  Indicator  with  LED 


pressure  displ^  and  a  0  -  5volt  analog  output.  An  Omega  Brand  K-lype  (Chromel- 
Alumel,  3.15  -  1445  K)  dimnocouple  connected  to  an  Omega  Brand  Model  DP25-TC-A 
programmable  digital  thennocouple  meter  disp%ed  the  tunnel  diamber  total  temperature 
and  provided  a  0  -  10  volt  analog  output,  with  a  scale  factor  of  283.15K/volt.  The 
uncertainties  in  die  total  pressure  and  temperature  were  ±  0.02  MPa  and  ±  2  K, 
respectively. 


Fltntr  Strtif|thlncr  ^  1 ; 
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Thick  Staiiilcvic  Steel  Itiictt 


(a)  Photographs 


(b)  Tunnel  Calibration  with  Old  Straightner 


(c)  Tunnel  Calibration  with  New  Straightner 
Figure  3.4  New  Flow  Straightner 


The  test  section  had  a  7.62  x  7.62  cn?  cross-section,  with  a  laigth  of  26.7  cm. 
The  test  section  sidewalls  were  constructed  from  type  347  stainless  steel  and  had  optical 
grade  glass  windows  (flat  to  within  52  nm/cm)  for  optical  measurements.  One  of  the  top 
walls  had  three  slots  with  corresponding  blocks  for  probe  insertion  at  different  axial 
positions,  and  the  other  was  a  solid  plate  manufactured  from  4.45  cm  thick  optical  grade 
Plexiglass  for  optical  access.  The  wind  tunnel  models  were  located  in  the  lower  wall. 
The  leading  edge  of  die  diamond  injectors  was  located  7.14  cm  downstream  of  the  nozzle 
exit  Shown  in  Fig.  3.5  is  die  test  section  geomeby. 


Fig.  3.5  Single-Port  Test  Section  Geometry  and  Coordinate  System 

The  variable  area  diffuser  was  7.62cm  wide  and  91  cm  long.  The  throat-area  was 
adjusted  to  7.6  cm  (28.3  cm  downstream  of  the  diffuser  inlet),  for  optimal  diffuser 
efficiency.  The  minimum  pressure  ratio  was  measured  at  21,  which  corresponds  to  a 
diffuser  efficien<y  of  approximately  75%.  Considering  the  compact  size  of  the  diffuser, 
this  efficiency  was  deemed  acceptable.  The  tunnel  air  discharged  to  die  atmosphere 
through  an  Allied  Witan  Co.  model  M60  air  sound  muffler,  which  was  attached  to  the 
diffuser  exit  through  a  15.24  cm  (6  inch)  pipe. 

S.l.S  Low  Reynolds  Number  Calibration  Wind  Tunnel 

The  cross- film  probes  are  sensitive  to  R^nolds  number.  Hence,  the  probes  had 
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to  be  calibrated  over  the  range  expected  across  the  injector  plumes.  Thus,  a  low 
Reynolds  number  calibration  in-draft  facility  was  constructed.  The  facility  consisted  of  a 
Mach  1.7  contoured  nozzle  connected  to  the  shock  tunnel  vacuum  receiver  tank.  The 
cross-section  of  the  calibration  facility  was  7.62  x  7.62  cn?.  Shown  in  Fig.  3.6  is  a 
photograph  of  the  calibration  facility.  The  measured  Mach  number  is  1.73,  and  die 
R^olds  number  per  meter  is  14  x  10*. 


Fig.  3.6  Low  Reynolds  Number  Calibration  Facility 


3.2  Wind  Tunnel  Models 

3.2.1  Single-port-injectors 

Four  diamond  injectors,  with  different  incidence  angles,  and  one  circular  injector 
(90°)  were  designed,  constmcted  and  tested.  Figure  3.7  shows  the  injector  schematics. 
The  effective  diameter  of  all  injectors  was  4.89  mm  The  ducts  for  die  diamond-shaped 
injectors  were  machined  to  angles  of  10°,  27.5°,  45°,  and  90°  respectively  with  respect  to 
the  surface  of  die  injector  exit  port.  A  135°  injection  was  achieved  by  operating  die  45° 
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injector  backwards.  All  of  the  single  port  diamond  injectors  had  a  15°  vertex  half  angle. 
This  small  angle  was  chosai  to  ensure  attached  shock  waves  near  the  wall  and  to 
minimize  the  effects  by  the  reflected  waves  from  the  side-walls.  The  injector  was 
mounted  on  die  bottom  plate  of  the  test  section  as  shown  in  Figure  3.5.  The  injectant  air 
was  induced  into  die  wind  tunnel  model  through  an  aluminum  duct  with  a  diameter  of 
1.27  cm  The  jet  exit  Mach  number  was  1.0.  The  jet-to-free  stream  momoitum  flux  ratio 
(J)  was  varied  by  varying  the  injector  total  pressure.  The  jet  flow  conditions  are  listed  in 
Table  3.2. 
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Fig.  3.8  Schematic  of  the  Single  Port  Injector  (45°  Injector  Shown) 
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Table  3.2  Test  Matrix 


Test 

Condition 

Injector 

geometry^ 

(deg) 

Cd‘ 

r 

V 

Pt 

(MPa) 

p/pJ^ 

1 

SD 

10.0 

0.23 

0.43 

0.08 

0.10 

8.4 

2 

SD 

27.5 

0.41 

0.43 

0.16 

0.10 

9.7 

3 

SD 

45.0 

0.67 

0.43 

0.27 

0.10 

10.0 

4 

SD 

90.0 

0.94 

0.43 

0.40 

0.10 

10.3 

5 

SC 

90.0 

0.99 

0.43 

0.42 

0.10 

10.3 

6 

SD 

135.0 

0.67 

0.43 

027 

0.10 

10.0 

7 

SDTR 

45.0 

0.67 

0.43 

027 

0.10 

10.0 

8 

SDDR 

45.0 

0.67 

0.43 

0.27 

0.10 

10.0 

9 

SD 

10.0 

0.24 

1.91 

0.41 

0.10 

41.1 

10 

SD 

27.5 

0.41 

1.91 

0.76 

0.45 

45.5 

11 

SD 

45.0 

0.68 

1.91 

1.29 

0.45 

46.4 

12 

SD 

90.0 

0.91 

1.91 

1.74 

0.45 

46.7 

13 

SC 

90.0 

0.90 

1.91 

1.72 

0.45 

47.3 

14 

SD 

135.0 

0.68 

1.91 

1.29 

0.45 

46.4 

15 

SD 

10.0 

0.15 

2.59 

0.33 

0.59 

54.6 

16 

SD 

27.5 

0.42 

2.59 

1.03 

0.59 

60.5 

17 

SC 

45.0 

0.67 

2.59 

1.68 

0.59 

61.7 

18 

SD 

90.0 

0.95 

2.59 

2.40 

0.59 

62.1 

19 

SC 

135.0 

0.67 

2.59 

1.68 

0.59 

61.7 

20 

SDTR 

45.0 

0.67 

2.59 

1.68 

0.59 

61.7 

21 

SDDR 

45.0 

0.67 

2.59 

1.68 

0.59 

61,7 

“SD:  single  port  diamond,  SC:  single  circular,  SDTR:  single  port  diamond  with  tapered  ramp, 
SDDR:  single  port  diamond  with  double  ramp.'’Injector  incidraice  angle.  -  ideal  momentum  ratio. 
Vo  =  effective  momentum  ratio  accounting  for  viscous  losses  [see  Bowersox  et  al  (2003)]. 
^Effective  pressure  ratio,  Tj/T..  =  4.0  and  m/«„  =  0.41 . 


The  tapered  ramp  of  Wilson,  Bowersox  and  Glawe  (1998)  was  chosen  to  enhance 
die  plume  vorticity,  l^ilere  in  the  Wilson  et  al.  stiufy,  this  ramp  produced  the  largest  gain 
in  plume  vorticity.  The  45°  diamond  injector  was  chosen  as  the  test  bed  for  the  variable 
wall  geometry  ramp.  This  ramp  was  designed  sudi  that  the  ramp  angle  matched  that  of 
the  injector  (45°)  and  die  hei^t  was  chosen  so  that  the  vertex  of  die  ramp  was  adjacent  to 
the  lower  side  of  the  plume.  These  measures  were  takai  to  ensure  strong  ramp-plume 
interaction,  ^)^hile  not  overwhelmingly  distorting  die  flow  widi  the  ramp.  Figure  3.8 
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shows  Ihe  final  ramp  geometry. 


Top  viCTV 


1.45  cm— ►! 


Front  view  Sidg>igw 


Fig.  3.9  Schematic  of  Variable  Wall  Tapered  Ramp 
[Design  Based  on  Wilson,  Bowersox  and  Glawe  (1998)]. 


3.3  Computational  Facilities 


The  bulk  of  the  numerical  simulations  were  performed  on  the  Air  Force  ASC 
HPC  Compaq  and  Origin  computer  systems. 


3.4  Instrumentation 

3.4.1  Surface  Oil  Flow  Visualization 

Dow  Coming  Brand  200  fluid  (Dimethicon)  with  a  kinematic  viscosity  of  1000 
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centistokes,  irtWch  is  a  medium  viscosity  linear  polydimettiylsiloxane  polymer,  was  used 
for  siuface  oil  flow  visualization.  Day-Glo  Color  T-15  blaze  orange  fluorescent  pigment 
was  mixed  with  the  fluid.  The  mixture  was  spread  on  the  bottom  surface  of  the  test 
section  and  illuminated  by  an  ultraviolet  lamp  through  the  top  surface  of  die  wind  tunnel, 
which  was  constructed  from  optical  grade  Plexiglass.  Photographs  of  the  oil  patterns 
were  taken  widi  a  Kodak  Professional  DCS  460c  color  digital  camera  that  was  mounted 
on  a  tripod  above  die  top  plate.  The  Kodak  camera  back,  which  was  attached  to  a  Nikon 
N90  camera  bo<fy,  had  a  high-resolution  3060x2036  pixel  color  Charge  Coupled  Device 
(CCD)  imager.  A  Nikon  N90  50  mm  lens  was  used,  and  the  shutter  speed  and  f  stop 
were  1.3  and  2.8,  respectively.  The  camera  was  connected  by  a  SCSI  port  and  cable  to  a 
NEC  personal  computer  with  a  300  MHz  Intel  Pentium  2  processor,  64MB  Sync  DRAM, 
8.4GB  hard  drive  PC  that  operated  with  Microsoft  Windows  98.  On  the  computer,  an 
Adaptec  AVA-1502  SCSI  Host  Adapter  with  a  SCSI2  port  transferred  the  digital  images 
recorded  by  CCD  imager  to  the  computer,  where  the  images  were  acquired  and  displayed 
through  the  Kodak  TWAIN  driver  into  the  Adobe  Photoshop  4.0  application  software  for 
Windows.  The  image  recording  was  controlled  via  the  Adobe  Photoshop  program. 

3. 4. 2  Two-Color  Digital  Particle  Image  Velocimetry  (PIV) 

The  laser  sources  used  for  the  present  PIV  measurements  were  two  Continuum 
Brand  Surelite  SL  I- 10  Nd:  NAG  laser  systems.  The  two  Surelite  Nd:NAG  lasers  each 
produced  coherent  green  (532  nm)  laser  beams.  The  red  (607  nm)  laser  beam  was 
formed  by  passing  one  green  beam  dirough  a  SulfaRhodamine  (640  nm)  methyl  alcohol 
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solution  dye  laser.  The  two  (green  and  red)  laser  beams  were  then  merged  in  space 
through  a  dichroic  splitter,  and  the  resulting  co- located  beams  were  passed  through  sheet¬ 
forming  optics  and  directed  into  the  test  section.  The  temporal  delay  between  the  red  and 
green  laser  was  set  at  200  ns  with  a  Stanford  Research  Systems  Model  DG535  four- 
channel  digital  delay/pulse  generator.  The  Kodak  Model  DCS  460  high  resolution  CCD 
camera  (described  in  the  previous  Section)  with  Nikon  N90  50mm  and  150  mm  lenses 
were  used  to  record  the  Mie  scattering  color  images. 

A  TSI  Brand  Model  9306  Six- Jet  Atomizer  was  used  to  produce  olive  oil  seed 
particles  with  mean  aerodynamic  diameter  <  0.6pm.  The  particles  were  injected  with  the 
injector  air.  Because  of  die  higji  pressures  associated  with  the  Mach  5.0  wind  tunnel,  the 
seed  generator  was  placed  inside  a  pressure  vessel. 

Because  of  the  hi^  fidelity  temporal  laser  separation  control,  the  digital  two- 
color  system  was  well  suited  for  hypersonic  flow.  The  PIV  system  provided  direct 
measurements  of  the  instantaneous  planar  velocity  (in  the  x  and  y  directions)  and 
vorticity.  In  addition,  the  instantaneous  Mie  scatter  images  provided  qualitative  insight 
into  die  large-scale  turbulent  structure  of  the  flow.  Unfortunately,  this  PIV  system  was 
not  robust  enough  to  acquire  enough  images  for  statisdcal  turbulence  measurements. 
More  specifically,  the  camera  acquisition  rate  translated  into  one  PIV  image  per  tunnel 
run.  Hence,  only  instantaneous  information  was  obtained. 

3.4.3  Mie  Scattering  Flew  Visualization 

Two  types  of  Mie  scattering  flow  visualization  images  were  acquired, 
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instantaneous  and  time  averaged.  Bar  the  instantaneous  images,  two  pulsed  Continuum 
brand  Surelite  Nd;YAG  were  used.  The  pulse  duration  was  10  ns  and  the  repetition  rate 
was  10  Hz.  The  lasers  produced  two  coherent  green  (532  nm)  beams.  One  beam  was  used 
to  pump  a  dye  [SulfaRhodamine  (640nm)  methyl  alcohol  solution]  laser  to  generate  a  red 
(640  nm)  beam.  The  two  (green  and  red)  laser  beams  were  merged  in  space  with  a 
dichroic  splitter,  and  the  resulting  co-located  beams  were  passed  through  sheet- forming 
optics  and  directed  into  the  test  section  through  an  optical  grade  Plexiglas  ceiling.  The 
temporal  delay  between  the  red  and  grerai  laser  was  set  at  200  ns  with  a  Stanford 
Research  Systems  Model  DG535  four-channel  digital  delay/pulse  generator.  The  two- 
color  system  allowed  for  instantaneous  velocity  measurement.  15  A  TSI  Brand  Model 
9306  six-jet  atomizer  was  used  to  seed  the  injectant  with  olive  oil  particles.  The  3dB 
frequency  response  of  for  this  0.5-  1.0  mircometer  diameter  seed  material  has  been 
previously  established  at  60  -  220  kHz.  19  The  scattered  light  from  both  the  red  and  green 
laser  pulses  was  collected  with  the  CCD  camera  described  in  the  previous  section.  For  the 
time  averaged  Mie  scattering  images,  only  one  laser  was  used  and  the  camera  exposure 
was  increased  so  that  30  laser  pulses  were  captured  on  a  single  CCD  image. 

3.4. 4  Pressure  Sensitive  Paint  (PSP) 

Pressure  sensitive  paint  (PSP)  provides  surface  pressme  distributions  with  larger 
areas  and  finer  spatial  resolution  than  conventional  pressure  ports.  The  ISS,  Inc.  Brand 
platinum  octaethylporphyrin  (PtOEP)  pressure  sensitive  paint  was  used  to  measure  the 
surface  pressure  distribution  around  the  jet  port.  It  had  absorption  maxima  at  excitation 
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wavelengths  of  380  and  540  nm,  with  an  emission  at  650  nm.  Because  the  present 
pressure  sensitive  paint  utilized  fluorescence,  a  relatively  high  frequency  response  (~  1 
msec)  was  achieved.  The  intensity  of  the  imaged  surface  is  governed  by  the  Stem- 
Volmer  relation  (i.e.,  /<//  =  i  +  cp/po),  where  c  is  constant  that  depends  on  the  molecular 
properties  of  the  paint  The  actual  calibration  was  slightly  more  general  (i.e.,  /«//  =  a  + 
bp/po). 

The  paint,  when  excited  with  a  blue  LED  light  source,  produced 
photoluminescence  in  order  to  lose  the  absorbed  energy,  and  some  of  the 
photoluminescence  was  quenched  by  the  local  oxygen  (a  pressure  dependent  parameter). 
Hence,  the  luminescence  intensity  was  pressure  dependent  The  luminescence  intensity 
generated  Ity  the  coating  was  then  collected  by  the  Pixelvision  SpectraVideo  Brand  CCD 
back-illuminated  grey-scale  camera,  model  SV512V1A/PFT-95  wifli  a  Nikon  frl6-50 
mm  lens.  Uncertainty  estimates  are  listed  in  Table  3.1. 

3.4.5  Traverse  and  Linear  Voltage  Displacement  Transducer 

The  traverse  system  consisted  of  a  15W  Bodine  Electric  Model  NSH-12  DC 
fractional  horsepower  gear  motor  that  drove  a  slide,  vsiiich  moved  the  probe  up  and  down 
in  the  j^direction.  Hie  traverse  speed  was  controlled  by  a  Minarik  Electric  Model  SL15 
controller.  The  speed  was  set  as  7.62  mm/s  (at  70%  of  full  power).  The  start  and  stop  of 
the  traverse  was  controlled  by  a  Magnecraft  Electric  W211ACPSOX-5  delay  on  start 
relay  and  a  Magnecraft  Electric  W21  lACPSRX-5  delay  on  ending  relay,  respectively.  To 
obtain  measurements  during  the  steady  flow  period,  the  traverse  system  was  deltyed  two 
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seconds  after  the  tunnel  start,  and  stopped  to  end  traverse  measurement  after  five 
seconds. 

The  position  of  the  probe  was  recorded  by  an  Omega  Model  LD 100- 150  Linear 
Voltage  Displacement  Transducer  (LVDT)  widi  nominal  linear  range  152.4  mm.  The 
LVDT  is  a  precision  variable  inductor  with  frequaicy  response  up  to  15  kHz.  An  Omega 
Model  SP200A  signal  processor  received  signals  from  LVDT,  and  generated  a  DC 
voltage  proportional  to  the  inductance,  which  was  linear  with  core  displacement.  The 
uncertainty  of  LVDT  measurement  was  0.067  mm. 

3.4. 6  Five-Hole  Pressure  Probe 

An  Aeroprobe  Inc.  Brand  five-hole-probe  was  used  to  measure  the  local  mean 
Mach  number,  Pitot  pressure,  static  pressure  and  total  pressure.  The  probe  was  a  19° 
semi  vertex  angle  axisymmetric  cone  constructed  from  stainless  steel,  with  a  blunted  tip 
containing  the  Pitot  probe.  There  were  four  elliptical  shaped  pressure  taps,  formed  from 
0.79  mm  outer  diameter  tube,  at  a  90°  interval  around  die  circumference  of  the  cone.  The 
cone  tip  was  soldiered  to  an  L-shaped  1.59mm  outer  diameter  stainless  tube,  which  held 
five  0.79  mm  outer  diameter  miniature  stainless  tubes.  The  space  between  the  tubes  was 
filled  with  steel  wire  and  epoxy  to  prevent  displacement.  The  L  shape  tube  was  inserted 
into  a  6.35mm  thick  and  14.6cm  long  stainless  streamlined  stmt  with  sharp  leading  and 
trailing  edges.  The  five  taps,  which  were  connected  to  the  five  0.79  mm  outer  diameter 
miniature  stainless  tubes,  were  connected  to  the  pressure  transducers  through  five  1.59 
mm  outer  diameter  miniature  Tygon  R-3603  tubes.  The  Pitot  pressure  was  measured  an 
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Omega  0-0.69  MPa  (0-100  psia)  Pxl39-100A4V  pressure  transducer,  and  the  four  static 
pressures  were  measured  with  four  Omega  0-0.10  MPa  (0-15  psia)  Pxl39-15A4V 
transducers.  All  five  transducers  were  excited  with  an  Omega  PSS-5A  power  supply. 
Each  transducer  had  an  output  0.25-4.25 volts  and  a  4.8mm  outer  diameter  inlet 
connector.  Based  on  the  normal  shock  and  Taylor-MacColI  conical  flow  theory,  the 
Mach  number  is  a  unique  function  of  the  ratio  of  the  cone-static  and  Pitot  pressures.  The 
numerical  conical  flow  and  Raleigh  Pitot  solution  is  tabulated  Table  3.2.  A  curve  fit  for 
the  Mach  number  as  a  function  of  the  pressure  ratio  is  given  by 


1  r  p 

--  =  -1.424  +24.873 
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V 


+1165.51 


The  percraitage  differences  between  fire  curve  fit  and  the  exact  theory  are  also  listed  in 
Table  3.2.  A  comparison  of  the  calibration  curve  fit  to  the  numerical  solution  is  given  in 
Fig.  3.11.  The  remaining  flow  properties  were  computed  using  the  usual  compressible 
flow  relationships  [Anderson  (1990)].  An  important  inference  from  Fig.  3.11  is  that  at 
high  Mach  numbers,  greater  than  5.0,  the  change  in  Mach  number  becomes  very  sensitive 
to  small  changes  in  the  pressure  ratio;  this  is  because  of  hypersonic  similarity.  Hence,  the 
uncertainties  in  the  Mach  number  are  amplified  at  higher  Mach  numbers.  The 
uncertainty  estimates  for  the  Mach  number  and  static  pressure  acquired  with  this  probe 
listed  in  Table  3.1  are  given  as  a  function  of  Mach  number.  It  is  also  noteworthy  that  the 
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present  study  was  focused  on  the  plume  and  boundary  layer  regions  of  the  flow  (M  <  5), 
and  thus,  die  hypersonic  similarity  difficulty  was  avoided. 

3.4. 7  Shadowgraph 

The  optical  system  for  the  shadowgraph  consisted  of  a  mercury  vapor  light  source 
and  a  mirror  to  form  a  parallel  column  of  light  The  beam  was  passed  through  the  test 
section  and  onto  a  reflective  screen  to  form  the  shadowgraph  image.  The  shadowgraph 
images  were  digitized  with  the  Kodak  CCD  camera  mentioned  above. 

3.4.8  Data  Acquisition  Boards  and  Testpoint  Software 

The  analog  signals  fi'om  die  various  probes  (except  the  cross- film,  discussed  in 
the  next  paragraph)  were  converted  to  digital  signals  with  a  12-bit  resolution,  Keithley 
Model  DAS-1802AO  low- gain  data  acquisition  board.  This  system  could  acquire  data  at 
a  rate  up  to  333  ksamples/s,  with  a  maximum  A/D  throughput  of  312.5  ksamples/s  for 
operation  in  the  0-5volt  unipolar  mode.  All  channel  inputs  were  widiin  a  range  of  05 
volt  in  order  to  increase  total  system  throughput  and  reduce  noise  by  using  the  same  gain 
for  all  channels.  The  system  is  software-configurable  for  sixteati  single-ended  onboard 
channels  or  eight  differential  channels.  The  data  board,  with  the  STA-1800U  screw 
terminal  I/O  interface  accessory,  was  connected  to  an  IBM  Pentium(r)  OverDrive  56MB 
RAM  PC  operating  Windows  98.  The  Capital  Equipment,  Inc.  data  acquisition  control 
software  package  Testpoint  was  used  to  control  the  data  acquisition.  The  total  number  of 
samples  and  sampling  frequency  were  10,000  and  1  kHz,  respectively.  Consequently,  the 
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data-taking  period  was  10  seconds. 


Flow  Solver,  Numerical  Domain  and  Solution  Methodology 

GASP  (Aerosofl,  Inc.)  is  a  commercial  finite-volume,  structured,  multi-block, 
compressible  code  widi  the  ability  to  solve  both  unsteady  and  steady  RANS  equations.  In 
additirm,  GASP  can  also  solve  the  parabolized  Navier-Stokes,  the  Euler  equations,  and 
the  incompressible  Navier-Stokes  equations.  Time  accurate  solutions  can  also  be 
obtained  through  an  implicit  or  explicit  procedure. 

Tuihulaice  can  be  modeled  in  GASP  using  a  variety  of  equations.  The  set  of 
equations  include  the  algebraic  model  of  Baldwin- Lomax,  die  one-equation  model  of 
Spalart-Allmaras,  and  two-equation  models  such  k-e  and  k-OD.  Inviscid  fluxes  can 
be  diosen  firom  the  standard  Roe,  Van  Leer,  Roe  widi  Hartrai’s  correction, 
preconditioned  Roe  flux,  AUSM+,  full  flux  for  solving  parabolized  Navier-Stokes,  and 
the  artificial  compressibility  method  for  solving  incompressible  flows. 

GASP  version  4.0  includes  new  features  such  as  the  ability  to  model  complex 
geometries  using  Chimera/overset  grids.  The  code  can  be  executed  on  shared  and 
distributed  memory  machines.  This  is  accomplished  by  using  the  message  passing 
interface  (MPI)  libraries. 

Initial  gimnlatinns  of  the  flow  were  performed  on  a  Cartesian  grid  in  which  die 
injector  was  modeled  using  a  stair  case  arrangement  of  cells  along  the  tunnel  wall.  A 
constant  velocity  profile  was  specified  as  die  boundary  condition  for  the  injector.  This 
was  found  to  be  inadequate  since  losses  in  the  injector  port  were  neglected  and  also  the 
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penetration  was  found  to  be  higher.  Availability  of  GASP  version  4.0  along  witii  the 
construction  of  a  twenty  four-node  Beowulf  class  computer  made  it  feasible  to  increase 
the  grid  size  to  obtain  better  resolution  of  the  flow.  Also  the  injector  port  was  modeled  to 
include  the  injector  port  losses. 

Figure  3.10  shows  flte  current  grid  used  in  die  simulations.  The  grid  is  composed 
of  six  blocks  with  a  total  of  5.26  x  10^  cells.  Grid  dimensions  are  listed  in  Table  3.3.  Only 
a  quarter  of  the  physical  domain  was  modeled  since  the  flow  is  symmetric  about  the 
centerline  of  the  tunnel  and  the  region  of  interest  was  below  a  height  of  yid  =  8.0.  Grid 
lines  are  clustered  near  the  tunnel  floor  to  satisfy  /  <  1.  Grid  lines  are  also  clustered 
close  to  die  injector  exit  to  capture  important  flow  features.  The  grid  lines  normal  to  the 
tunnel  floor  are  aligned  in  the  shock  direction  within  the  limits  of  the  desired  grid  quality. 
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Figure  3.10  Computational  Grid 


Table  3.3  Grid  Dimensions 


Block 

Fine  Grid 

Medium  Grid 

Coarse  Grid 

1 

21x151x101 

11x76x51 

6  X  26  X  26 

2 

281  X  151  X  101 

141x76x51 

141x26x26 

3 

81  X  151  x41 

41x76x21 

41x26x11 

4 

21x151x31 

11x76x16 

11x26x16 

5 

81x51x41 

41x26x21 

41x26x11 

6 

21x51x31 

11x26x16 

11x26x16 

Total 

5260000 

657500 

118125 

Cells 

Inlet  conditions  for  the  tunnel  were  obtained  from  a  flat-plate  calculation  to  match 
e?q)erimental  boundary  layer  thickness  5.  Adiabatic  wall  conditions  were  used  for  die 
tunnel  floor  and  the  injector  wall.  Symmetry  conditions  were  used  for  the  plane  along  the 
centerline  of  the  tunnel.  Density  and  pressure  were  used  for  tfie  injector  inlet  conditions. 
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Slip  boundary  conditions  were  used  for  the  outer  wall  of  the  tunnel.  All  other  boundaries 
were  extrapolated.  Figure  3.11  shows  the  computational  domain  with  associated 
boundary  conditions. 


Figure  3.11.  Boundary  Conditions. 

The  k-o)  [Wilcox  (2000))  turbulence  model  with  compressibility  correction  was 
used  in  this  simulation.  Roe  flux  [Roe  (1981)]  along  with  the  min-mod  [HBrsch  (1992)] 
limiter  and  a  third-order  upwind  scheme  was  used  for  die  invisdd  flux,  laminar  transport 
properties  were  modeled  using  die  Sudierland  viscosity  and  conductivity  models  and  die 
constant  Schmidt  number  Difilisivity  model.  A  turbulent  Prandd  number  of  0.9  and 
turbulent  Schmidt  number  of  0.5  were  used.  Gauss-Seidel  technique  along  widi  Gauss- 
Seidel  inner  iterations  was  used  to  obtain  the  solution  for  all  global  iterations. 

Both  grid  and  solution  convergaice  were  based  on  the  caiterline  floor  pressures. 
Grid  conveigence  was  performed  in  situ  using  the  grid  sequencing  technique  available  in 
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The  number  of  cells  in  eEdi  gnd  level  is  listed  in  Table  3.3.  The  centerline 
pressures  for  the  three  grid  levels  are  given  in  Fig.  3.12.  Solution  conveigoice  was 
monitored  by  plotting  the  centeriine  pressures  every  25  iterations.  The  solution  was 
deemed  converged  when  the  change  in  die  centerline  pressure  was  negligible.  Figure  3.12 
shows  a  representative  plot  of  the  centerline  pressure  for  the  90°,  0.1  MPa  jet  total 
pressure. 


Figure  3.12  Solution  Convergence  Plot 
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Chapter  4 

Results  and  Discussion 


4.1  Experimental  Results  and  Discussion 

Surface  Oil  Flow  Visualization 

Surface  oil  flow  visualizations  were  acquired  for  all  of  the  no-ramp  conditions 
listed  in  Table  3.2.  The  results  are  shown  in  Fig.  4.1.  The  diamond  injector  results  were 
divided  into  two  basic  categories:  attached  and  detached  leading  edge  shocks. 
Representative  examples  of  each  diamond  injector  type  are  shown  in  Figs.  4.2a  and  4.2b, 
where  the  freestream  is  from  left  to  right,  and  the  jet  flow  is  directed  outward  from  the 
paper.  Shown  in  Fig.  4.2c  is  the  high-pressure  circular  injector  result.  Examining  the 
attached  case  given  in  Fig.  4.2a  (low-pressure,  a  =  27.5°),  the  first  flow  feature,  moving 
from  left  to  right  is  a  curved  line  of  oil  accumulation  enclosing  the  diamond-shaped 
injector  exit.  This  line  is  attached  to  the  injector  leading  edge,  and  is  essentially  co¬ 
located  with  the  jet  interaction  shock  in  the  region  near  the  injector;  this  is  confirmed 
with  the  pressure  sensitive  paint  measurements  described  below  in  the  Pressure  Sensitive 
Paint  section.  However,  further  downstream,  the  shock  wave  strength  decreased  and 
approached  a  Mach  wave,  and  hence  it  was  no  longer  coincident  with  die  line  of 
separation.  The  second  curved  line  of  oil  accumulation  that  wrapped  around  the  jet 
indicates  the  location  of  a  small  horseshoe  vortex.  Third,  two  wake  vortex  lines  due  to 
H)-wash*®  were  observed  to  emanate  at  the  injector  trailing  edge.  Fourth,  between  the 
wake  vortex  lines  is  die  reattachment  line.  In  general,  the  flow  topologies  for  the 
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detached  shock  cases  (e.g..  Fig.  4.2b)  were  similar  to  those  of  the  attached  shock  case. 
The  obvious  exceptions  are  (1)  die  detached  leading  edge  jet  interaction  shock  and  (2)  the 
horseshoe  vortex  was  more  pronounced  and  it  reached  furflier  upstream  for  die  detached 
shock  case.  The  circular  injector  flow  topolo©?  (Fig.  4.2c)  was  very  similar  to  the 
detached  diamond  injector  (Fig.  4.2b).  However,  when  compared  to  the  detached 
diamond  injector  flow,  die  circular  injector  shock  standoff  distance  was  larger  and  the 
angle  betweai  the  wake  vortex  lines  was  lower.  In  general,  the  topology  described  above 
is  very  similar  to  the  lower  Mach  number  results  of  Tomioka  et  al.*® 


Fig.  4.1  Summary  of  Surface  Oil  Flow  Visualizations 
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(a)  Low-Pressure,  a  =27.5° 


(c)  High-Pressure,  a  =90° 

Fig.  4.2  Representative  Surface  Oil  Flow  Visualizations 


Upon  examination  of  the  surface  oil  flow  visualizations  for  all  of  the  injection 
angles  (10  -  90°)  at  both  injection  pressures,  the  following  specific  observations  were 
drawn.  First,  the  state  of  the  leading  edge  shock  wave  (attadied  or  detached)  for  tiie 
diamond  injectors  dep^ded  on  both  a  and  the  injector  pressure.  The  shocks  w^ 
attached  for  the  low-pressure,  a  =  10°  and  27.5°  cases.  The  shock  standoff  distance  was 
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very  small  for  the  45°  case.  As  a  increased  from  10°  to  45°,  the  leading  edge  shock  angle 
increased  from  35°  to  52°.  The  increase  is  explained  by  recognizing  that  the  discharge 
coefficient,  and  hence  mass  flux,  increased  significantly  with  injector  angle  (see  Table 
3.2).  Moreover,  as  a  increased,  the  3-D  relieving  effect^®  decreased. 

Focusing  now  on  die  medium-pressure  results,  the  shocks  were  attached  for  a  = 
10°.  For  a  =  10°  to  45°,  the  shock  angles  were  nominally  10°  larger  dian  die 
corresponding  low-pressure  cases.  Also,  compared  to  the  low-pressure  data,  the  high- 
pressure  wake  regions  were  wider  in  the  z-direction;  the  wakes  and  wake  vortices 
extended  further  downstream,  and  the  horseshoe  vortices  reached  further  upstream.  For 
the  detached  shock  cases,  the  high-pressure  shock  strength  was  larger  than  die 
corresponding  low-pressure  results,  as  indicated  by  the  larger  lateral  extent  of  the  nearly 
normal  shock  region  and  die  larger  shock  wave  detachment  distance.  The  explanation  for 
all  of  these  observations  is  that  the  hi^er-pressure  jets  expanded  more  wdien  injected  into 
the  freestream  creating  additional  blockage  and  penetration,  and  thus  the  above  listed 
observations.  Similar  observations  hold  for  the  highest  pressure. 

Shadowgraph  Photography 

Shadowgraph  images  were  acquired  for  all  of  the  no-ramp  conditions  listed  in 
Table  3.2.  The  results  are  summarized  in  Fig.  4.3.  In  these  figures,  the  cross  flow  is  from 
left  to  right,  and  the  injector  port  is  in  die  tunnel  floor.  The  present  discussion  will  focus 
on  the  three  cases  described  above.  First,  the  attached  shock  case  corresponds  to  the  low- 
pressure,  a  =  27.5°  test  condition.  A  clearly  defined  jet  interaction  shock  is  depicted,  and 
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the  barrel  shocks  are  faintly  visible.  Whai  compared  to  the  medium  pressure,  a  =  90 
case,  the  main  difference  between  die  detached  and  attached  shock  shadowgr^hs  is  the 
presence  of  the  lambda  shock,  \nhich  was  created  by  the  stronger  detached  bow  shock 
wave  interaction  with  the  tunnel  boundary  I^ct.  The  shock  structure  for  the  circular 
injector  was  similar  to  that  for  the  detached  injector. 


Fig.  43  Summary  of  Shadowgraph  Images 


Examination  of  all  of  die  iiadowgraphs  confirmed  the  trends  described  above. 
Specifically  for  the  low-pressure,  a  =  10°  and  27.5°  cases,  leading  edge  lambda  shocks 
were  not  present  indicating  an  absence  of  upstream  flow  separation.  For  a  =  10°  to  45°, 
the  shock  angles  increased  from  30°  to  49°  with  jet  incidence  angle,  demonstrating  diat 
the  shock  strength,  and  hence  total  pressure  loss,  increased  with  jet  incidence  angle. 
Quantification  of  the  shock  losses  is  described  in  the  later.  The  90°  diamond  and  90° 


46 


circular  cases  both  had  lambda  shocks,  indicating  flow  separation.  Barrel  shocks  were 
visible  for  the  45°  -  90°  cases.  However,  th^  were  not  visible  for  the  remaining  two  lower 
angle  cases.  It  was  expected  that  for  the  lower  angles,  the  barrel  shocks  were  not  as 
strong  due  to  the  lower  discharge  coefficients,  and  subsequent  reduced  pressure  (listed  in 
Table  3.2).  In  addition,  the  jets  did  not  penetrate  veiy  far  into  die  freestream;  hence,  the 
barrel  shocks  were  blurred  the  tunnel  boundary  layer.  The  initial  shock  angles  for  each 
incidence  angle  were  noticeably  larger  (8°  -  13°)  for  the  high-pressure  conditions  as 
compared  to  the  corresponding  low-pressure  cases.  Lambda  shocks  were  present  for  all 
but  the  lowest  angle  case,  and  the  barrel  shocks  were  clearly  visible  in  all  of  the  high- 
pressure  shadowgraphs. 

Mie  Scattering  Imagery 

Summarized  in  Fig.  4.4  are  the  mean  Mie-scattering  flow  visualizations. 
Compared  in  Figs.  4.5a  and  b  are  the  instantaneous  images  for  the  low-pressure,  a  =  27.5 
and  90®  cases.  In  both  images,  the  freestream  flow  is  fi*om  left  to  right.  The  instantaneous 
image  shows  that  large-scale  turbulent  structures  were  present  in  the  jet  plumes;  however, 
the  size  of  the  stmctures  near  the  jet  orifice  appeared  to  have  increased  with  increasing 
incidence.  For  the  90®  case  shown  in  Fig.  4.5b,  the  tunnel  was  intaitionally  run 
approximately  5.0K  below  the  oxygen  liquefaction  temperature  to  produce  Rayleigh 
scattering  fi-om  the  condaised  liquid  OTQ^gen  in  the  freestream.  Hence,  the  tunnel 
boundary  layer,  just  upstream  of  the  injector  port  and  die  jet  interaction  shockweare 
visible  due  to  vaporization  of  the  oxygen  in  the  warmer  boundary  layer  and  across  the 
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shock.  This  rather  unique  flow  visualization  indicated  fliat  evai  fliough  the  jet  interaction 
shock  for  tins  case  was  detadied  near  fbe  wall  (Table  2),  it  re-attached  to  flie  jet  near  the 
point  where  the  plume  bait  over  towards  the  fireestream  direction,  hi  addition,  the  jet 
plume  narrowed  in  this  region,  which  is  indicative  of  the  plume  spreading  more  in  die 
lateral  direction;*^  Fig.  5a  shows  a  similar  phenomenon. 


Low-Pnanrt  HlgM>re«Mire 


fig.  4.4  Mean  Mie  Scattering  Image  Summary 


The  time-averaged  data  (wcanqiles  are  shown  in  Figs.  4.5c-e)  show  more  clearly,  the 
upper  edge  and  haice  mean  penetration  of  the  plume  into  the  fi'eestream  Time  averaged 
Mie  scattering  images  were  obtained  for  the  first  two  pressure  cases  listed  in  Table  3.2.  It 
was  also  apparent  that  the  plume  narrowing  described  in  the  previous  paragraph  was 
present  for  all  of  the  diamond  injectors,  and  absent  for  die  circular  injector.  Visual 
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(a)  Upper:  Low-Pressure,  a  =  45° 


(b)  Lower:  Low-Pressure,  a  =  90° 


(c)  Low-Pressure,  a  =27.5° 


h-Pressure.  a  =  90° 


(e)  High-Pressure,  a  =  90° 

Fig.  4.5  Representative  Mie-Scattering  Flow  Visualizations 
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measurements  of  the  centerline  penetration  fyp/d)  were  obtained  from  these  data  for  x/d 
less  than  about  8.0.  The  jagged  plume  edges  seen  in  Fig.  4.5  resulted  in  scatter  in  the  yp/d 
measurements  ranging  from  about  ±0.1,  close  to  the  jet  exits,  to  approximately  ±0.2  at  the 

downstream  limit  of  die  measurement  (x/d  «  8.0). 

Various  power-law  correlation  methods  have  been  proposed  for  penetration.  ’  A 
generalized  form  for  these  correlations  for  the  present  sonic  injection  into  a  Mach  5.0 
freestream  was  writtrai  as  yp/d  =  $(§)»  where  ^  =  Jeff''’(x/dy .  For  the  present  study,  Jeff  ~ 
Jsina,  and  the  values  of  J  are  listed  in  Table  2.  Plotted  in  Fig.  4.6  are  all  of  the 
penetration  data  with  the  above  scaling.  As  indicated  in  the  plot,  the  diamond  and  circular 
injector  data  collapse  reasonably  well  (variation  in  yp/d  is  within  the  scatter  of  the  data) 
onto  to  two  separate  curves,  one  for  the  diamond  injectors  and  the  other  for  the  circular 
injectors.  For  both  data  sets,  n  was  found  to  be  nominally  0.33;  this  value  is  consistent 
with  previous  studies  at  lower  supersonic  freestream  Mach  numbers.^^  Values  of  0.15  and 
0.26  for  m  were  found  to  collapse  the  diamond  and  circular  port  data,  respectively.  The 
function  i^(^)  for  die  circular  injector  data  was  fte  straight  line  yp/d  =  1.93|±0.40,  wfrere 
the  linear  form  agrees  widi  previous  studies.  The  nonzero  intercept  reflects  the  fact  that 
the  plume  Kcpanded  upstream  of  the  injector  port  (see  Fig.  4.5d).  For  the  diamond 
injector,  ^(^)  was  linear  to  witiiin  the  experimental  scatter  for  ^  greater  than  about  1.0, 
with  a  slope  of  2.65.  However,  for  ^  less  than  1.0,  the  data  systematically  dipped  below 
the  linear  trend.  A  hyperbolic  curve  was  incorporated  to  fit  the  data  over  die  full  range  of 
^  shown  in  Fig.  6;  the  resulting  equation  was  yp/d  =  [H-(2.65^)i^]^^-l.  The  reproduction 
of  the  penetration  was  within  the  expected  uncertainty  of  die  data 
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Penetration  measurements  at  x/d  =  7.5  are  listed  in  Table  2  (column  8)  for  all  of 
the  test  conditions.  From  these  data,  die  following  specific  conclusions  were  drawn.  First, 
for  the  diamond  injectors,  the  penetration  at  x/d  =  7.5  increased  by  70  and  90%  as  a 
increased  from  10  to  90°,  for  the  low-  and  high-pressure  cases,  respectively.  Second,  for 
given  incidence,  the  penetration  increased  by  nominally  30%  going  from  low  to  hi^ 
pressure  for  all  of  the  diamond  injectors;  for  die  circular  injector,  the  increase  was  50%. 
Also,  comparing  the  90°  diamond  and  circular  injectors,  the  magnitude  of  the  penetration 
was  essentially  the  same  at  x/d  =  7.5.  The  data  and  correlations  shown  in  Fig.  4.6  also 
show  that  the  penetration  for  the  circular  injectors  at  both  pressures  exceeded  that  of  the 
90°  diamond  injectors,  for  x/d  <  7.5.  However,  the  plume  trajectories,  i.e.,  slopes,  for  die 
90°  diamond  injectors  were  approximately  30%  larger  than  those  for  the  circular  port  at 
x/d  =  7.5.  Hence,  it  is  anticipated  that  the  downstream  penetration  for  the  90°  diamond 
injectors  would  exceed  that  of  the  circular  injectors,  which  agrees  with  a  previous 
study.** 
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Fig.  4.6  Penetration  Correlation  Results 
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Mach  Number  Contour  Data 


Fig  4.7  Mach  number  Contour  Summary 
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The  Mach  number  contours  are  shown  in  Fig.  4.7.  Examining  the  Mach  number 
contour  for  the  low-pressure,  a  =  45°  diamond  shaped  injector  first  (enlarged  contour  on 
die  left).  Starting  at  the  top,  the  interaction  shock  wave  is  the  roughly  diamond  shaped 
structure  with  the  vertex  located  at  z/d  «  0.0,  y/d  «  4.2,  and  that  ardied  down  to  z/d  » 
±3.9  at  y/d  «  2.0.  The  narrow  high-Mach  region  below  the  shock  is  the  invisdd  shock 
layer  between  the  shock  and  the  plume.  The  lower  extent  of  the  shock  layer  corresponds 
to  the  upper  boundary  of  the  plume,  which  at  the  centerline  (z/d  =  0.0)  was  located  at  y/d 
«  2.8-3.2.  Now  focusing  on  the  bottom  of  the  contour,  the  wake  or  up-wastf  ®  region,  is 
seen  as  the  low  Mach,  semicircular  region  centered  at  z/d  =  0.0.  The  vertical  extent  of 
this  region  (y/d  «  1.4)  coincides  with  lower  “edge”  of  the  plume  as  inferred  from  the  Mie 
scattering  images.  The  tunnel  boundary  layers  are  also  noticeable  on  both  sides  of  the 
contour.  The  boundary  layer  thinning  in  the  areas  adjacait  to  the  wake  resulted  from 
entrainment  of  low  momentum  boundary  layer  fluid  into  the  plume  and  wake.^*  Above 
the  wake  region  is  the  plume.  Aldiougji  the  plume  edges  are  not  discemable  in  Fig.  4.7, 
the  expected  kidney  bean®  shape  of  the  plume  can  be  inferred  from  the  Mach  contours. 
Specifically,  starting  at  the  midpoint  of  the  lower  edge  of  the  plume  (z/d  =  0.0,  y/d  « 
1.4),  and  traversing  around  the  plume  in  a  counterclockwise  direction,  the  plume 
boundary  is  expected  to  have  wr^ped  down  and  around  the  wake  region,  dipping  closer 
to  the  wall  in  the  region  where  the  boundary  layer  thinned  (z/d  »  1.0).^*  Although,  the 
exact  lateral  extent  of  the  plume  is  not  discernible,  the  “dog-ear”  structure  located  near 
z/d  »  1.0,  y/d  «  1.5  is  indicative  of  the  low- momentum  boundary  layer  fluid  induction 
into  the  right  vortex  associated  with  the  counter-rotating  vortex  pair.^®’*^  Symmetry  about 


53 


the  z-axis,  to  wifliin  the  experimental  xmcertainty,  was  observed. 

The  Mach  number  contour  for  the  detached-shock  diamond  injectors  are  also  in 
Figs.  4.7.  The  overall  interpretation  of  these  images  is  similar  to  that  of  the  attached 
shocks.  The  notable  differences  are  (1)  the  interaction  shock  extended  about  further  into 
the  freestream;  (2)  the  shock  shapes  had  more  curvature,  which  is  consistent  with  a 
detached  bow  shock,  and  (3)  the  low-Mach  upwash  region  was  significantly  larger 
indicative  of  die  increased  mass  flow  rate  and  penetration 

The  Mach  number  contour  for  the  circular  injector  is  given  in  Figs.  4.7,  where  the 
features  are  very  similar  to  the  detached  diamond  injectors.  However,  the  penetration  of 
the  interaction  shock  was  approximately  10%  lower  than  that  for  corresponding  90- 
degree  diamond  injector.  Also,  the  circular  injector  low-Mach  upwash  region  was 
significantly  smaller  than  that  for  the  diamond  injector.  However,  the  dog-ear 
structures  at  z/d  *  i  1.0,  y/d  ^  1.5  are  more  pronounced  for  the  circular  injector,  which 
indicates  stronger  plume  vorticily. 

Pressure  Sensitive  Paint 

Pressure  Sensitive  Paint  (PSP)  data  were  acquired  for  all  of  the  no-ramp  conditions 
listed  in  Table  3.2;  the  results  are  summarized  in  Fig.  4.8.  The  pressure  distributions 
reflect  the  flow  topologies  described  with  surface  oil  flows  described  earlier.  Focusing  on 
the  lowpressure,  45°  injector  as  a  representative  attached  shock  case,  Ihe  pressure  ratio 
(p/pbo)  increased  to  approximately  2.3  across  die  interaction  shock.  The  surface  pressure 
subsequently  decreased  through  an  expansion  as  the  flow  wrapped  around  the  injected 
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fluid  to  a  relatively  low  wake  pressure,  p/p®  »  0.3;  this  low  pressiu’e  was  expected.*  For  a 
represaitative  detached  diamond  injector  case,  consider  the  higji-pressure,  90°  injector 
surface  pressure  distribution.  For  diis  case,  the  horseshow  vortex,  the  high-pressure 
structure  that  wr^ped  around  the  injector,  was  a  dominant  high'prcssure  flow  structure. 
The  peak  pressure  ratio  within  the  horseshoe  vortex  was  ^proximately  4.2.  The  pressure 
ratio  in  the  wake,  just  downstream  of  the  injection  port  was  nominally  0.3.  The  surface 
pressure  field  for  the  circular  injector  contained  the  same  structures  as  for  die 
corresponding  90-degree  diamond  injector.  However,  die  shock  standoff  distance  was 
significantly  larger  (123%)  than  for  the  corresponding  diamond  flow.  The  peak  pressure 
ratio  for  the  hi^pressure  circular  injector  was  4.9;  die  minimum  pressure  ratio  in  die 
wake  was  0.15.  The  effective  backpressure  was  measured  fi-om  the  surface  pressure  plots 
as  die  average  pressure  around  die  periphery  (within  ^proximately  0.2d  of  the  injector 
edges)  of  the  injector,  die  resulting  ratios  (A.)  fi'om  this  analysis  are  listed  in  Table  4.1. 


Table  4.1  Performance  Summary 


a* 

f) 

m 

II 

Pj^a 

17 

10.0 

1.9 

1.2 

7.0 

0.99 

27.5 

2.6 

1.3 

7.5 

0,95 

45.0 

3.0 

1.6 

6.3 

0.89 

90.0 

3.4 

2.0 

5.4 

0.78 

90.0" 

2.0 

5.4 

0.71 

10.0 

2.5 

1.4 

29.4 

0.94 

27.5 

3.5 

1.7 

26.8 

0.84 

45.0 

4.1 

1.7 

27.3 

0.74 

90.0 

4.7 

2.0 

23.4 

0.70 

90.0** 

4.8 

2.2 

0.69 

“First  five  entries  correspond  to  P,j  -  0.  lOMPa  (±1 .0%);  the  second 
five  correspond  ioPf  =  0.46  MPa  (±3.0%),  ’’Circular  Injector 
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For  perpendicular  injection  iiroiigh  circular  orifices,  the  initial  jet  penetration  has 
been  shown  to  be  strongly  depend^t  on  tiie  R/peb-^"^  Haice,  that  pressure  ratio  was  also 
computed  and  included  in  Table  4.1.  As  indicated,  for  each  of  die  two  constant  total 
pressure  data  sets,  pj/peb  decreased  by  25-30%  as  a  increased  fi-om  10  to  90®.  It  is 
expected^'*  that  this  decrease  worked  to  dcCTease  the  initial  penetration,  and  thus 
counteracted  a  portion  of  tiie  penetration  increase  that  would  have  resulted  from  tiie 
increase  in  Jeff  that  resulted  fi:om  the  increase  in  incidence  angle. 

Centeriine  pressure  traces  for  all  of  the  cases  in  Table  2  were  obtained  fiom  the  PSP 
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measurements  and  are  plotted  in  Fig.  4.9.  TTie  leading  edge  of  the  injector  was  at  x/d  = 
0.0.  The  trailing  edge  of  the  diamond  and  circular  injectors  were  at  x/d  =  2.4  and  1.0, 
respectively.  The  peak  pressure  just  upstream  of  the  injector  for  the  detached  cases  and  a 
given  geometiy,  increased  with  increasing  injector  pressure.  For  a  given  pressure  ratio, 
the  90°  circular  injectors  produced  larger  peak  pressures  than  did  the  90°  diamond 
injector,  indicating  stronger  shocks  for  the  circular  injectors.  Standoff  distances  for  the 
detached  shock  waves  as  estimated  from  this  figure  are  summarized  in  Table  2.  The 
standoff  distance  increased  with  both  injector  pressure  and  incidence  angle,  and  was 
largest  for  the  circular  injector.  The  wake  pressure  traces,  from  just  downstream  of  the 
injector  ports  to  x/d  «  3.5,  for  all  but  the  low-pressure  10°  diamond  injector,  were  similar 
in  magnitude  and  shqie.  The  low-pressure,  10°,  diamond  injector  had  a  veiy  low 
discharge  coefficient,  and  hence  the  overall  flow  disturbances  were  small.  The  gap  in  the 
data  between  x/d  =  3.5  to  5.0  was  the  result  of  a  imperfection  in  the  paint  coverage  due 
to  the  sealing  of  a  port  in  the  floor;  this  port  is  clearly  visible  in  the  contours  shown  in 
Fig.  8.  Downstream  of  this  region,  the  flow  was  compressed,  and  subsequently  relaxed 
back  to  nearly  the  freestream  pressure. 
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Fig.  4.9  Centerline  Pressure  Traces 


Approximate  Shock  Loss  Analysis 

In  order  to  estimate  the  bow  shock  total  pressure  loss,  the  three-dimensional 
interaction  shock  shape  was  modeled  from  the  surface  shock  shtqie,  tiie  center  plane 
shock  shape,  and  cross-sectional  shock  shape.  Hie  cross-sectional  shock  shape  was 
estimated  from  tiie  five-hole  probe,  low-pressure  Madi  number  data  at  x/d  =  8.0.  The 
cross  sectional  jet  interaction  shock  shapes,  nondimensionalized  by  the  maximum  y,  z 
coordinates,  at  x/d  =  8.0,  were  fit  with  the  following  3”*  order  polynomial 

=«(^/^max)^  +  +  «^(^/ ’  \ihere  tiio  coofficients  a-d  were  found 

by  applying  boundary  conditions  at  z/^nax  =  0  and  z/%iax  =  1.  The  best  agreemrait  to  the 
data  was  achieved  with  the  following  boundary  conditions:  At  z/^ax  =  0,  y^max  =  1  and 
d(y/ymax  )/d(z/4,ax)  =  -0.4.  At  z/%iax  =  1,  y/ymax  =  0  and  d(y/ymax  )/d(z/^aax)  =  -0.3.  The 
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mean  square  errors  for  all  cases  varied  from  4.0  -  7.0%.  The  shock  shape  along  the  tuimel 
surface  [i.e.,  ^ax(x/d)]  was  evaluated  from  the  data  in  Figs.  4.3  and  4.8.  The  caiterline 
shock  shq)e  [i.e.,  ymax(x/d)]  was  directly  measured  from  the  shadowgr^h  photographs 
(Fig.  4).  With  die  3-D  shock  sh^e  modeled,  the  unit  normal  vector  (n )  to  the  shock 
surface  was  calculated  by  computing  the  gradient  to  die  shock  shape  function.  With  the 
unit  normal,  the  normal  component  of  the  Mach  number  over  die  shock  surface  was 
found  by  the  relation  M„=Mjl -n.  Because  die  jet  was  introduced  into  the  tunnel 

boundary  l^er,  Mi  was  evaluated  as  a  function  y/d  using  die  Van  Driest  n  conqiressible 
turbulent  velocity  profile  scaling.^^  The  shock  total  pressure  was  then  estimated  with  the 
normal  component  of  the  Mach  number  using  normal  shock  dieory.^"  The  resulting  area 
averaged  total  pressure  losses,  that  is 


s 


(4.1) 

were  computed  for  x/d  <  8,  where  S  is  the  shock  surface  area,  and  die  integrand  in  die 
numerator  is  a  function  of  Mn  for  air.  The  results  are  listed  in  Table  2.  For  the  higji- 
pressure  cases,  the  normalized  shock  cross-section  fit  was  assumed.  In  order  to  quantify 
the  scatter  due  to  die  curve  fit,  die  shock  loss  was  also  calculated  for  the  45°  case  using 
two  additional  curves  that  defined  upper  and  lower  bounds  in  the  experimental  data  for 
the  scaled  cross  sectional  shock  fit.  The  resulting  uncertainty  was  approximately  ±4.0%. 
The  shock  loss  data  in  Table  4.1  demonstrate  diat  losses  decreased  widi  decreasing  jet 
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incidence  angle  and  injection  pressure.  The  largest  pressure  losses  were  incurred  by  Ae 
90°,  circular  injector. 

4.2  Numerical  Results  (GASP  v  4.1) 

The  numerical  results  obtained  wiA  Ae  GASP  code,  as  described  m  Chapter  3, 
are  Ascussed  m  this  section.  The  salimt  parameters  described  above  provide  Ae  frame 
work  for  Ae  initial  Ascussion.  However,  Ae  last  two  topics  describe  in^rtant  results 
that  do  not  have  experimental  counterpart  data. 

Separation  Distance 

The  second  column  of  Table  4.2  lists  Ae  computed  separation  Astances 
normalized  by  Ae  effective  Aameter  of  Ae  jet.  The  results  mAcate  Aat  Ae  separation 
Astance  increases  wiA  mcreasmg  /  and  mcreasmg  angle  of  mcideace  Of.  The  27.5°,  low 
momentum  case  produces  an  attached  shock  while  all  oAer  shocks  are  detached. 
Experimental  results  show  Aat  m  boA  Ae  27.5°  and  45°,  low  momentum  case  and  27.5°, 
hi^  momentum  case  Ae  shocks  are  attached.  Figure  4.10  shows  Ae  representative 
numerical  oil  smear  images  for  all  Ae  momentum  ratios,  for  Ae  45°  mjection  case.  The 
freestream  flow  is  from  rigfit  to  left  m  Ae  figure.  The  jet  mteractioo  shock  is  shown 
clearly  for  all  three  momentum  ratio  cases.  It  is  at  Ae  location  ^ere  Ae  freestream  flAd 
deviates  from  Ae  main  flow  Arection.  The  shock  standoff  Astance  mcreases  wiA 
mcrease  m  Ae  momentum  ratio.  The  leading  edge  horseshoe  vortex  angle  also  mcreases 
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with  increase  in  the  mommtum  ratio.  This  leads  to  an  increase  in  the  reattachment 
distance  as  indicated  in  the  figure.  Also,  tiiere  is  an  increase  in  die  lateral  distance 
betweai  the  tunnel  centerline  and  the  outer  boundary  of  the  leading  edge  horseshoe 
vortex.  This  increase  is  also  noted  in  die  Lateral  Penetration  section.  The  separation 
between  the  wake  vortices  also  increases  with  increasing  momentum  ratio.  The  trend  of 
increasing  shock  separation  distance  with  increase  in  momentum  ratio  is  observed  in  all 
cases. 


Table  1.2  Numerical  Simulation  Results. 


a* 

(deg) 

y/d’ 

y/<f 
%  Error 

n% 

Error 

27.5 

AS 

2.10 

0.93 

19.23 

2.11 

45 

0.19 

3.00 

0.88 

0.00 

1.12 

90 

0.70 

2.00 

0.80 

41.18 

2.56 

135 

2.00 

1.85 

0.85 

NA 

NA 

27.5 

0.26 

4.24 

0.80 

21.14 

4.76 

45 

0.46 

5.36 

0.68 

30.73 

8.11 

90 

1.90 

3,30 

0.47 

19J9 

32.86 

135 

6.80 

3,22 

0.50 

NA 

NA 

27.5 

0.30 

4.75 

0.76 

NA 

NA 

45 

0.58 

5.63 

0.62 

NA 

NA 

90 

3.60 

4.10 

0.39 

NA 

NA 

135 

9.75 

4.08 

0.44 

NA 

NA 

“First  four  entries  correspond  toP^  =  O.lOMPa;  the  second  four  correspond  toP^  -  0.46  MPa;  the  third 
four  correspond  toP,y  =  0.63  MPa,  **AS  =  Attached  Shock,  =  7.5,  NA=Not  Applicable 
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c.  Ptj  =  0.63  MPa 

Fig.  4.10  45°  Injection -Numerical  Oil  Smear. 


Normal  Penetration 

The  experimental  penetration  data  ei,  x  Id  =  7.5  listed  in  Table  4. 1  indicates  that 
the  penetration  increased  widi  increasing  angle  of  inddaice  and  increasing  momentum 
ratios.  Normal  poietration  of  the  injectant  fluid  in  the  current  simulation  was  measured 


by  projecting  the  injectant  streamlines  on  tiie  symmetry  plane  of  the  tunnel.  Figure  4.11 
shows  the  injectant  streamlines  for  the  45°  cases.  The  cross  plane  is  at  a  location 
ofxld=1.5.  Both  the  cross  plane  and  the  streamlines  are  colored  with  Mach  number 
values.  The  penetration  height  was  taken  at  the  streamline  with  the  highest  y  Id  value. 
The  data  summarized  in  Table  4.2  indicated  that  the  penetration  distance  increases  by 
almost  80  percent  when  the  total  pressure  of  the  jet  increases  from  0.1  MPa  to  0.46  MPa. 
However  the  increase  in  jet  total  pressure  from  0.46  MPa  to  0.63  MPa  results  in  only  a  5 
percent  increase  in  paietration.  A  similar  trend  is  observed  in  the  27.5°  case. 

Lateral  Penetration 

Figure  4.12  shows  the  lateral  penetration  of  the  injectant  fluid  streamlines.  It  is 
observed  that  although  the  penetration  in  the  normal  direction  seems  considerably  higher 
for  the  high  jet  total  pressure  cases,  it  is  not  as  uniform  as  the  low  jet  total  pressure  case. 
The  lateral  penetration  shows  a  marked  increase  as  the  jet  total  pressure  is  increased  and 
the  penetration  is  more  uniform  vriien  compared  with  the  normal  penetration.  Also,  clear 
increase  in  the  normal  and  lateral  penetration  is  noted  as  flie  jet  total  pressure  is  further 
increased  to  0.63  MPa. 


63 


c  =  0.63  MPa 

Fig.  4.11  45**  Injection  -Injectant  Streamlines. 


(a)  P^i  =  0.1  MPa  (b)  P,j  =  0.46  MPa  (c)  P^  *  0.63  MPa 

Fig.  4.12  45“  Injection  — Injectant  Streamlines  Shosring  Lateral  Penetration. 

Figures  4.13-4.15  show  the  injectant  streamlines  for  tfie  27.5®,  90®,  and  135® 
cases.  The  core  of  the  jet  fluid  penetrates  less  laterally  in  the  27.5®  cases  as  compared  to 
the  45®  cases.  It  is  seen  that  the  27.5°  injection  follows  a  trend  similar  to  flie  45°  cases 
with  increase  in  jet  total  pressure.  The  90®  and  135®  injection  cases  show  marked  increase 
in  lateral  penetration,  especially  at  higher  momentum  ratios.  The  injectant  travels  far 
upstream  in  the  135®,  high  momaitum  ratio  case  as  seen  in  Fig.  4.15b.  This  is  not 
observed  in  other  cases. 

Total  Pressure  Loss 

Total  pressure  loss  parameter  [Eq.  (4.1)]  was  calculated  for  all  the  cases.  A  plane 
at  station  xld  =  8.0,  spanning  ztd  =  -4  to  +4  and  y  Id  =  0  to  7  was  used  for  evaluating 
the  total  pressure  after  the  shock.  The  resulting  values  are  listed  in  Table  4.2.  It  is  seen 
that  as  the  angle  of  incidence  increases  from  27.5®  to  90®  die  total  pressure  loss  increases. 
This  closely  parallels  the  experimental  results  in  Table  4.1.  This  increase  in  total  pressure 
loss  can  be  attributed  to  stronger  shocks  that  are  produced  as  the  injector  angle  increases. 
The  differraice  between  the  experimental  and  numerical  total  pressure  loss  increases  for 
cases  where  the  jet  total  pressure  is  0.46  Mpa 


The  135°  injection  cases’  total  pressure  loss  is  less  than  tiiat  of  tiie  90°  cases.  This 
is  contrary  to  the  expected  result  since  a  135°  injection  should  result  in  a  stronger  shock. 
Additional  investigation  into  this  phaiomenon  needs  to  be  performed. 


a  Pm  =  0.1  MPa 


b.  P,j  =  0.46  MPa 

Fig.  4.13  27.5*  Injection  -Injectant  Streamlines. 
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b.Pq  =  0.46  MPa 

Fig.  4.15  135”  Injection  -Injectant  Streamlines. 


Vortex  Structures 

The  transverse  jet  interaction  flow  field  is  dominated  1^  numerous  vortices  as 
described  in  Chapter  2  In  addition  to  previously  documented  vortices,  a  new  vortex 
mechanism  is  observed  in  the  current  simulations.  Figure  4.16  shows  the  vortex  cores 
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[Haimes  and  Kaiwright  (1999),  Sujuki  and  Haimes  (1995)]  in  Ae  flow  field.  Starting 
from  Ae  bottom  left  comer  of  Ae  figure  and  moving  to  Ae  top  ri^t  comer  we  first 
encounter  Ae  vortex  core  for  Ae  shock  separation  vortex.  Above  Ae  trailing  edge  of  Ae 
jet  are  Ae  vortex  cores  for  Ae  newly  identified  transverse  counter  rotating  vortex  pair 
(TCVP).  Starting  near  Ae  TCVP  are  Ae  vortex  cores  for  Ae  wake  vortices.  The 
reattachmoit  vortec  cores  start  near  Ae  trailing  edge  of  Ae  jet. 

Figure  4.17  shows  some  of  Ae  vortices  associated  wiA  Ae  cores  Ascussed 
previously.  Hie  leading  edge  shock  separation  vortex  is  clearly  visible  at  Ae  bottom  left 
comer  of  Ae  figure.  Moving  along  a  Aagonal  paA  we  raicounter  Ae  TCVP.  This  vortex 
pair  is  formed  as  Ae  mcoming  boundary  layer  fluid,  which  is  transported  mto  Ae  upper 
l^ers,  interacts  wiA  Ae  fluid  issuing  from  Ae  jet  ©dt 


Fig  4.17  Transverse  Counter  Rotating  Vortex  Pair. 
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Chapter  5 

Conclusions 


5.1  Experimental  Conclusions 

An  experimental  study  to  characterize  the  near- field  fy/d  <  8.0)  flow  for  gaseous 
injection  throu^  diamond  shaped  orifices  (15°  half  angle)  at  four  incidence  angles  (10, 
27.5,  45  and  90°)  and  three  total  pressures  (Py  =  0.10,  0.46  and  0.63  MPa)  into  a  high 
Reynolds  number  {Re/m  =  53  x  10®)  Madi  5.0  fi-eestream  was  performed.  A  90°  circular 
injector,  with  the  same  exit  port  area  and  total  pressures,  was  examined  for  comparative 
purposes.  Flow  control  ramps  [Wilson,  Bowersox  and  GHawe  (1998)]  were  investigated 
to  quantify  the  effect  on  plume  vorticily.  The  experimental  methods  included  surface  oil 
flow  visualization,  shadowgraph  photography,  Mie-scattering  flow  visualization,  pressure 
sensitive  paint  (PSP)  and  a  Pitot-caie  five-hole  pressure  probe.  Flowfield  docum^itadon, 
jet  poietration  and  shock  induced  total  pressure  loss  were  derived  from  these  data.  From 
diese  data,  the  following  specific  conclusions  were  drawn.  First,  near-wall  attadbment  of 
the  jet  interaction  shock  wave  depeided  on  both  incidaice  angle  and  injector  pressure, 
where  the  shock  remained  attached  at  higher  incidence  angles  at  the  lower  injection 
pressire.  Second,  5ld  =  1.6;  hence,  die  detached  shocks  separated  the  incoming 
boundary  layer  and  generated  die  familiar  lambda  shock  structure.  Lambda  shocks  were 
not  observed  for  the  attached  shock  cases.  Third,  near-field,  x/d  <  8.0,  centerline 
penetration  correlations  yp/d  =fp(^,  where  ^  =  J^(x/d)”,  were  developed  for  both  the 
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diamond  and  circular  injectors;  fp  was  hyperbolic  and  linear  for  die  diamond  and  circular 
injectors,  respectively.  Fourth,  for  the  detached  cases,  the  highest  surface  pressures  \^ere 
associated  with  the  horseshoe  vortex.  Fifth,  the  upstream  peak  pressures  for  die  detached 
cases,  for  a  given  injector  geometry,  increased  with  increasing  injector  pressure.  Sixth, 
for  a  given  pressure  ratio,  the  90°  circular  injectors  produced  largest  peak  pressures. 
Seventh,  the  shock  induced  total  pressure  losses  decreased  with  decreasing  jet  incidence 
angle  and  injection  pressure,  where  die  largest  losses  were  incurred  by  die  90°,  circular 
injector. 

5.2  Numerical  Conclusions 

Numerical  simulations  were  performed  for  sonic  injection  through  diamond 
orifices  into  a  Mach  5.0  ffeestream  Four  different  injector  angles  (27.5°,  45°,  90°  and 
135°)  and  diree  total  pressures  (0.1,  0.46  and  0.63  MPa)  resulting  in  a  total  of  12  cases 
were  executed.  The  derived  data  from  these  simulations  were  compared  to  experimental 
results.  Although  some  of  die  data  compared  well  with  experimental  results,  diere  were 
significant  differences  in  other  cases.  Shock  separation  distance  was  found  to  depend  on 
the  jet  incidence  ai^le  and  die  jet  total  pressure.  Increasing  angle  of  incidence  and 
increasing  total  pressure  led  to  larger  shock  standoff  distances.  Penetration  was  found  to 
be  highest  for  the  45°  injection  as  compared  to  the  90°  case  (experimental).  But  the  90° 
injection  case  penetration  was  found  to  be  more  uniform  as  compared  to  the  45°  case. 
Lateral  penetration  results  indicated  that  higher  jet  total  pressures  results  in  more  lateral 
penetration  of  the  injectant  fluid.  Injectant  in  the  high  momentum  ratio,  135°  case 


72 


reached  far  upstream  of  the  injector  location.  Total  pressure  loss  follows  experimental 
trends  for  the  27.5°,  45°  and  the  90°  cases  while  the  total  pressure  loss  for  the  135°  cases 
is  less  compared  to  die  90°  cases.  Vortex  core  plots  showed  the  existence  of  a  transverse 
counter  rotating  vortex  pair  immediately  downstream  of  the  jet. 

5  J  Ongoing  Calculations 

High-fidelity  numerical  simulations  are  continuing  at  Texas  A&M  University 
using  Mentor’s  SST  model  with  detached  eddy  simulation  The  focus  of  these 
simulations  is  improved  understanding  of  the  vorticity  dynamics  and  the  controllability  of 
the  secondary  flow  motion  using  gas  dynamic  actuation.  Computer  support  for  these 
simulations  is  being  provided  by  the  Air  Force  ASC  HPC.  Student  support  for  analyzing 
the  numerical  data  is  being  provided  by  Texas  A&M  University.  The  results  from  these 
simulations  will  be  reported  to  the  Air  Force  and  AFOSR  at  a  later  date. 
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